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I.  Introduction 

One  of  the  principal  prerequisites  in  the  development  of  high  current, 
charged  particle  accelerators  and  high  power  lasers  is  the  development  of  a 
high  voltage,  low  inductance  nanosecond  switch.  This  requirement  has  led  to 
the  investigation  of  breakdown  phenomena  in  high  pressure,  overvolted  gaps. 
When  a  high  voltage  is  applied  to  a  parallel  plate  gap,  containing  a  gas  or 
mixture  of  gases,  the  ionization  processes  resulting  from  the  creation  of 
some  initial  electrons  in  the  gap  space,  and  the  subsequent  motion  of  the 
charge  carriers,  can  cause  the  collapse  of  the  voltage  across  the  gap.  Thus, 
the  evolution  of  the  ionization  in  the  applied  electric  field,  from  the  small 
number  of  initiatory  electrons  up  to  a  final  steady  current,  is  of  fundamental 
inportance  in  the  development  of  fast,  high  power  gas  switches. 

This  research  program  is  responsible  for  investigating,  theoretically 
and  experimentally,  the  gas  breakdown  phenomenon  as  a  function  of  overvoltage 
(i.e.  voltages  above  D.  C.  breakdown  voltage)  and  to  develop  a  theoretical 
model  for  the  observed  phenomena. 

In  the  following  sections,  we  discuss  the  findings  of  this  study. 


II.  Theoretical  Investigations 


There  is  at  present  no  agreement  among  researchers  as  to  what  processes 
are  important  in  the  breakdown  phenomenon  as  the  applied  voltage  is  increased 
above  self breakdown.  As  part  of  this  program,  a  colloquium  was  given  at  the 
naval  surface  weapons  Center,  Dahlgren;  in  which  the  various  theories  of  break¬ 
down  and  the  present  status  of  the  subject  was  reviewed.  As  a  result,  a  review 
paper  was  published  and  is  included  in  this  report  as  Appendix  I. 

In  our  research  effort,  a  model  to  describe  the  initial  phase  of  breakdown 
in  high  pressure,  overvolted,  spark  gaps  was  developed.  The  model  offers  the 
attractive  feature  of  the  unification  of  all  the  proposed  breakdown  phenomena; 
that  is,  it  provides  a  continuous  picture  of  the  development  of  the  initial 
phase  of  breakdown  above  the  Townsend  regime.  A  paper  describing  this  model 
has  been  published.  A  copy  is  included  as  Appendix  II.  A  conference  paper 
was  also  given  and  a  copy  of  the  abstract  is  also  included. 

A  mathematical  formulation  of  the  above  model  is  in  progress.  The  math¬ 
ematical  problem  is  basically  that  of  solving  the  coupled  Boltzmann-Poisson 
system  of  equations  in  a  region  of  high  electric  field.  This  is  to  be  done 
computationally.  A  major  obstacle  in  the  simulations  is  the  computational 
time  required  to  solve  the  Poisson  equation  in  cylindrical  coordinates.  A 
fast  alogarithm  is  required  in  order  to  be  able  to  tract  the  development  of 
the  avalanche  for  a  long  period  of  time.  As  a  result,  we  developed  a  fast 
alogarithm  for  solving  Poisson's  equation.  The  outline  of  a  paper,  which  is 
in  preparation,  and  the  computer  printout  of  the  alogarithm  (which  is  a  sub¬ 
routine  in  the  main  program)  is  included  in  Appendix  III.  An  explanation  of 
how  the  subroutine  is  used  is  also  included. 


Ill  Experimental  Investigations 


A  photograph  of  the  facility  built  for  the  experimental  study  of  gas 
breakdown  at  high  overvoltages  is  shown  in  Fig.  1.  A  schematic  layout  of 
this  facility  is  shown  in  Fig.  2.  The  setup  consists  of  three  main  parts: 

(1)  The  pulse  generator  (a  modified  Heds  pulser^),  (2)  the  experimental 
chamber  and  (3)  the  diagnostic  probes. 

The  Heds  pulser  is  capable  of  producing  trapezoidal  voltage  pulses  in 
the  range  of  50  to  160  kV  amplitude,  with  a  4  nsec  risetime  and  a  50  nsec 
width,  into  a  52  ohm  line. 

The  experimental  chamber  is  a  coaxial  transmission  line  with  52  ohm  charact¬ 
eristic  impedance,  terminated  in  a  matched  load.  An  interruption  in  the 
center  conductor  of  the  line  is  used  as  the  spark  gap. 

The  salient  features  of  this  facility  are: 

a)  Well  matched  components  to  minimize  reflections 

b)  Base  chamber  pressure  in  the  10”^  mmHg  range  to  minimize  unwanted 
impurities. 

c)  ability  to  use  large  electrode  separation  (up  to  5  cm)  while  main¬ 
taining  strong  and  relatively  homogenous  external  electric  fields. 

d)  direct  access  to  either  cathode  or  anode  (depending  on  the  polarity 
of  the  incident  pulse,  which  can  be  reversed. 

e)  capacitive  and  B  probes  at  various  locations  to  measure  incident, 
reflected  and  transmitted  pulses. 

Using  the  set  up  described  above,  we  have  carried  out  an  investigation 
of  the  observational  delay  time,  i.e.  the  time  elapsed  from  the  application 
of  the  voltage  to  the  gap  to  the  time  of  rapid  current  growth.  We  have  obtained 
data  of  the  behavior  of  this  time  lag  as  a  function  of  percent  overvoltage 


Fig.  1  View  of  the  experimental  facility 


nsec  UV  flash 


Schemotic  of  Experimental  Facility. 


(ranging  from  10%  to  1300%  above  D.C.  breakdown),  gas  impurity  (Technical 
grade  nitrogen,  98%  N2,  and  purified  nitrogen,  99.998%  N2),  and  pressures 
(ranging  from  50  to  1500  torr).  The  experiments  were  carried  out  in  a  1  cm 
gap,  using  highly  polished,  Rogowskii  shaped.  Aluminum  electrodes.  The 
observational  time  lag  consists  of  the  statistical  time  lag  and  the  formative 
time  lag.  The  role  played  by  the  statistical  phenomena  in  these  experiments 
is  different  than  for  small  gap  experiments.  This  is  due  to  the  fact  that 
for  our  case,  the  inter-electrode  volume  (large  gap  separation)  is  large,  so 
that  there  is  a  higher  probability  of  having  an  electron  present  in  the  gap 
at  the  time  the  voltage  is  applied  to  the  gap. 

Since  the  percent  overvoltage  is  a  parameter  of  importance  in  our  studies, 

a  D.C.  breakdown  study  was  carried  out  to  determine  the  self  breakdown  voltage, 

Vg.  These  experiments  were  done  using  aluminum  electrodes  and  two  types  of 

gases:  a)  Technical  grade  Nitrogen  (99.0%  N2),  and  b)  "purified"  Nitrogen 

(99.998%  N2)-  Examples  of  vs.  pd  curves  are  shown  in  Fig.  3  where  we 

2 

also  plot  the  results  obtained  by  Dakin  et.al.  D.C.  Breakdown  was  "defined" 
to  be  the  voltage  at  which  consistent  D.C.  breakdown  was  observed.  This 
voltage  however  has  a  very  narrow  distribution  so  that  changes  in  the  "defi¬ 
nition"  yields  little  difference. 

In  the  pulse-breakdown  study,  for  a  given  incident  pulse,  the  reflected 
and  transmitted  pulses,  due  to  the  collapse  of  the  gap,  were  recorded.  From 
the  transmitted  pulse  recordings,  the  observational  time  lag  can  be  determined. 
A  typical  recording  is  shown  in  Fig.  4.  The  part  of  the  signal  indicated  with 
an  "A"  is  the  capacitive  coupled  signal,  and  corresponds  to  the  arrival  of 
the  incident  pulse  to  the  gap.  A  time  later,  the  current  rapidly  increases, 
indicative  of  the  complete  collapse  of  the  insulating  properties  of  the  gas. 
This  elapsed  time  is  the  observational  time  lag.  Since  the  slope  of  the 


Fig.  3  Static  Breakdown 
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Fig.  4  Transmitted  current  pulse 


current  rise  is  very  high,  the  observational  time  lag  can  be  measured  to 
within  .1  nanosecond  accuracy;  using  a  Tektronix  7104  oscilloscope.  An 
expanded  display  of  the  region  "A"  is  shown  in  Fig.  5.  Note  that  the  current 
rise  is  so  high  that  the  scope  cannot  trace  it.  The  tail  end  of  the  trans¬ 
mitted  pulse  can  be  seen  on  the  right  of  the  photograph. 

Approximately  90  shots  were  taken  at  each  operating  point  (i.e.  for  a 
given  pressure  and  %  overvoltage).  The  procedure  used  in  gathering  this  data 
is  as  follows:  the  chamber  is  initially  evacuated  to  approximately  10'®  Torr 
and  then  filled  with  N2  to  a  pressure  of  a  few  Torr.  At  this  pressure,  the 
gap  is  broken  down  15-20  times.  The  chamber  is  once  again  evacuated  and 
back  filled  to  the  desired  operating  pressure.  Although  15  shots  might  not 
be  enough  to  thoroughly  "clean"  the  electrodes,  the  sparks  observed  at  the 
higher  pressures  are  randomly  distributed  over  the  electrode  surface  indicating 
statistical  independence.  Moreover,  the  number  of  electrons  emitted  per  second 
were  determined  to  be  relatively  stable  throughout  the  experiment  (this  point 
will  be  discussed  later). 

From  the  90  shots  taken  per  point,  a  histogram  of  the  observational  delay 
is  compiled.  A  typical  histogram  is  shown  in  Fig.  6.  A  running  average  delay 
time  is  also  tabulated  and  plotted.  As  seen  from  Fig.  7,  the  running  average 
delay  time  reaches  its  asymptotic  value  after  approximately  35  shots.  Obvious¬ 
ly  if  the  histogram  is  narrower,  the  asymptote  is  reached  earlier. 

To  describe  the  breakdown  statistics,  let  Idt  be  the  probability  that 
an  electron  will  appear  in  the  gap,  in  the  time  interval  between  t  and  t  +  dt. 

I  is  the  probability  per  unit  time  that  electrons  are  emitted  from  the  cath¬ 
ode,  or  the  emission  rate  (assuming  all  other  processes  that  lead  to  electron 
generation  are  small.  This  point  is  yet  to  be  justified  in  our  case  since  we 
have  a  large  interelectrode  volume).  Moreover,  let  W  be  the  probability  that 
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Fig.  5  Capacitive  coupled  signal  observed  prior  to  the  colapse  of  the  gap 
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Fig.  6  Histogram  of  breakdown  delay  time 
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Fig.  7  Running  average  delay  time 
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once  an  electron  appears  in  the  gap,  it  will  lead  to  breakdown.  These  proba¬ 
bilities  are  assumed  to  be  independent  of  shot  number. 

Knowing  the  above  probabilities,  it  is  then  possible  to  calculate  the 
probability  that  breakdown  will  occur  with  delay  time  t  or  greater.  This 
can  be  done  by  noting  that: 

[Probability  that  breakdown  does  not  occur  in  time  t  +  dt]  =  [Probability 
that  breakdown  does  not  occur  in  time  t]  x  [Probability  it  does  not  occur 
in  time  dt]. 

If  Pg  is  this  probability,  then 

Pg(t  +  dt)  =  Pg(t)  [1  -  WIdt] 
expanding  and  integrating  yields: 

PB(t)  = 

Pg(t)  can  be  determined  experimentally  as  the  ratio  N/N^,  where  N  is  the 
number  of  breakdowns  observed  to  have  a  delay  t  or  greater  and  is  the 
total  number  of  breakdown. 

Thus, 


or 

lln  (N/Njj)l  =  Wit 

Therefore,  a  plot  of  the  above  equation  will  yield  a  straight  line,  whose  slope 
is  WI.  Such  plot  is  shown  in  Fig.  8.  From  this  figure  we  determine  WI 

p 

to  be  -  5  X  10°  electrons  per  second.  If  we  assume  that  W  =  1,  then  the 

o 

cathode  emission  rate  is  5  x  10  electrons/second. 

Theoretically,  the  intersection  of  this  straight  line  with  the  t  axis 
gives  the  minimum  delay  time;  i.e.  the  formative  time  lag.  From  Fig.  8, 
we  note  that  there  is  a  deviation  of  the  experimental  curve  from  a  straight 
line.  This  deviation  is  a  function  of  the  width  of  the  distribution.  In 


Fig.  9,  a  broader  distribution  is  shown.  The  corresponding  Probability 
plot  is  shown  in  Fig.  10.  Note  that  the  deviation  is  greater.  However, 

O 

the  slope  of  the  straight  section  is  2.5  x  10  electrons/second  which  is 
approximately  equal  to  that  in  Fig.  8.  Further  studies  are  necessary  to 
complete  the  picture.  Moreover,  these  results  should  be  compared  to  small 
gap  experiments,  where  I  is  closely  associated  with  cathode  emission,  and 
to  UV  illuminated  experiments  where  I  is  very  large. 

A  plot  of  the  average  delay  time  as  a  function  of  percent  overvoltage 
is  shown  in  Fig.  11.  The  minimum  observational  time  lag  seems  to  also  follow 
the  same  shape.  However,  further  investigations  are  necessary  to  confirm 
this  behavior. 

The  fact  that  the  average  delay  time  has  a  minimum  as  a  function  of 
percent  overvoltage  was  rather  puzzling  at  first.  However,  it  can  be  quali¬ 
tatively  understood  from  our  model.  We  can  show  that  there  is  an  optimum 
value  of  percent  overvoltage  for  which  space  charge  distortion  of  the  field 
in  the  gap  is  maximum  thus  allowing  for  the  production  of  "runaway"  (from 
the  avalanche  head)  electrons  which  we  have  postulated  to  be  the  fundamental 
mechanism  for  avalanche  propagation.  A  quantitative  analysis  needs  to  be 
carried  outT  however,  further  data  is  necessary  (i.e.  UV  illumination  and 
small  gap  results  to  get  a  better  understanding  of  I  and  W). 

These  results  were  presented  at  the  33rd  Gaseous  Electronics  Conference. 
The  abstract  of  the  paper  is  included  in  Appendix  IV. 
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Fig.  11  Average  delay  time  vs  E/N 
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Electrical  Breakdown  of  Gases:  The  Prebreakdown 

Stage 

E.  E.  KUNHARDT,  member,  ieee 


Abstract-Xxi  this  paper,  a  review  of  the  theories  and  experiments 
devoted  to  the  understanding  of  the  development  of  the  electrical 
breakdown  of  a  gas  insulated  gap,  i.e.,  the  switching  delay,  is  presented. 
The  presentation  is  chronological.  The  classical  Townsend  and  streamer 
models  for  breakdown  are  discussed;  followed  by  a  brief  account  of  the 
continuous  acceleration  and  avalanche*  chain  models.  These  last  two 
models  have  been  proposed  primarily  to  describe  breakdown  at  large 
electric  Helds,  Then,  the  two-group  model  for  breakdown  at  volt¬ 
ages  above  approximately  20-percent  self-breakdown  is  presented. 
Finally,  a  brief  analysis  is  given  of  the  present  state  of  the  field  and  the 
direction  it  is  taking. 

I.  Introduction 

HE  PHYSICS  OF  HIGH  energy  densities  has  recently 
received  a  great  deal  of  attention  as  a  common  and  crucial 
area  of  interest  for  scientists  working  on  high-power  lasers, 
fusion,  high-current  charged-particle  accelerators,  and  weapons- 
effect  simulators.  In  the  present  scenario,  the  research  pro¬ 
grams  in  these  areas  rely  heavily  on  pulse -power  technology  to 
meet  their  goals.  One  of  the  principal  problems  limiting  the 
further  growth  of  this  technology  is  the  development  of  a 
switching  device  that  will  allow  the  fast  and  repetitive  transfer 
of  energy  from  an  energy  storage  device  to  various  transducers 
(e.g.,  laser,  particle-beam  generator,  etc.).  To  meet  these 
diverse  requirements  [1]  a  number  of  novel  switching  con¬ 
cepts  have  been  proposed.  In  many  of  these  approaches, 
switching  is  accomplished  by  causing  a  normally  insulating  gas 
to  undergo  a  transition  to  a  conducting  state.  The  various 
devices  that  operate  in  this  fashion  differ  mainly  in  the  way 
this  transition  is  initiated  and  in  the  characteristics  of  the  final 
conducting  stage,  i.e.,  whether  it  be  a  glow  or  arc  discharge. 
The  two  most  popular  devices  that  belong  to  this  category  of 
switches  are  the  thyratron  tube  and  the  pressurized  spark  gap. 
Presently,  these  devices  are  capable  of  operating  in  the  range 
of  kilovolts  and  kiloamperes. 

In  these  switching  devices,  transition  from  the  off  to  the  on 
state  is,  in  practice,  thought  of  consisting  of  two  phases:  1)  the 
switch  delay,  i.e.,  the  time  elapsed  from  the  desired  initiation 
of  the  transition  to  the  10-percent  value  of  the  current;  and  2) 
the  current  rise  time,  i.e.,  the  10-percent  to  90-percent  peak 
current  time.  Correlation  between  these  two  time  intervals 
and  the  physical  processes  occurring  in  the  switch  depend  on 
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the  particular  switching  device.  However,  one  can  say  tha 
general,  the  delay  time  is  determined  by  the  prebreakdc 
stage  of  the  insulating  gas,  whereas  the  current  rise  timi 
determined  by  a  number  of  factors:  the  dynamics  of  the  f 
conducting  state,  the  geometry  of  the  switch,  and  possi 
external  circuit  constraints.  In  this  paper,  we  would  like 
present  a  review  of  the  work  related  to  the  prebreakdc 
stage.  The  purpose  is  not  to  compile  a  catalog  of  referent 
but  to  give  an  overview  of  the  field  and  its  trend  which 
novice  can  use  as  a  starting  point  for  further  research.  1 
exposition  of  the  material  will  be  presented  chronologica 
This  will  help  the  newcomer  appreciate  the  present  state 
understanding,  and  review  for  the  expert  the  critical  steps 
the  development  of  the  field. 

Much  of  the  recent  research  effort  has  gone  into  minimiz 
the  current  rise  time,  and  we  are  now  approaching  a  lii 
where  the  total  transition  time  is  due  mostly  to  the  delay  tii 
The  reduction  of  the  switch  delay  time  will  require  a  bet 
understanding  of  the  physics  of  the  prebreakdown  stage.  T 
physics  is,  to  a  degree,  common  to  all  gas  switches,  and  ‘ 
been  very  difficult  to  investigate.  This  difficulty  stems  fro 
first,  the  broad  parameter  space  traversed  by  the  gas  as 
undergoes  breakdown;  second,  the  large  number  of  inter 
tions  taking  place  in  the  gas;  and  third,  the  fact  that  all  fi 
variables  are  inhomogeneous  in  space.  This  is  further  cc 
pounded  by  the  time  scales  involved  in  this  progression,  i 
in  the  nanosecond  range  in  some  cases.  This  requirement 
speed  and  the  sensitivity  constraints  imposed  on  the  eqi 
ment  have,  until  very  recently,  prevented  the  experimer 
investigation  of  many  of  the  prebreakdown  processors.  Ne 
less  to  say,  investigations  have  concentrated  on  idealized  (fr 
a  practical  point  of  view)  situations. 

As  an  illustration  of  these  difficulties,  the  ionized  com 
nent  of  the  gas  first  behaves  as  individual  charged  partit 
interacting  with  a  background  of  neutrals,  then  as  a  nonic 
plasma  (i.e.,  Debye  length  r£)>  plasma  dimensions)  : 
finally  as  an  ideal  plasma  (i.e.,  Tq  <  plasma  dimensions), 
single  analytical  framework  can  be  used  to  model  this  tra 
tion.  Moreover,  lack  of  light  emission  from  the  very  e; 
stages  and  the  short  time  scale  involved  has  prevented  meas’ 
ment  of  the  space-time  evolution  of  the  electron  number  c 
sity  during  this  stage. 

The  text  is  divided  into  several  sections.  In  the  first  sec^ 
a  standard  terminology  is  established.  The  subsequent  sect: 
form  part  of  one  continuous  story;  however,  they  have  I: 
separated  so  as  to  single  out  major  events  in  the  history  of 
field.  I  conclude  with  some  comments  on  the  present  ti 
in  the  field. 


0093-381 3/80/0900-0 130$00.75  ©  1980  IEEE 


KUNHARDT:  ELECTRICAL  BREAKDOWN  OF  GASES 


25 


11.  Statement  of  the  Problem 

The  physical  problem  to  be  discussed  is  primarily  the  growth 
of  the  ionization  of  a  gas  in  an  electric  field  and  the  subse¬ 
quent  breakdown  of  the  insulating  properties  of  the  gas  (a  very 
comprehensive  and  up-to-date  source  of  information  on  this 
subject  is  given  in  [2] ). 

Two  experimental  approaches  have  been  used  to  study  this 
phenomenon.  In  the  first  approach  [3] -[5]  a  dc-electric  field 
is  applied  between  two  electrodes  in  a  gas  (see  Fig.  1).  The 
cathode  is  irradiated  with  a  beam  of  ultraviolet  radiation,  gen¬ 
erating  a  constant  photoelectric  current  Iq  at  the  cathode.  For 
constant  values  of  £/p,  where  E  is  the  applied  electric  field  and 
p  is  the  pressure  of  the  gas  between  the  electrodes,  measure¬ 
ments  of  the  field  intensified  current  /  as  a  function  of  elec¬ 
trode  separation  d  are  made  (Fig.  2).  For  a  given  i’/p,  a 
separation  is  reached  at  which  the  current  becomes  self- 
sustaining.  That  is,  the  current  I  will  continue  to  flow  in  the 
circuit  even  after  the  external  source  of  ionization  has  been 
terminated.  This  self-sustaining  state  is  called  a  Townsend 
discharge  [6] .  Under  these  conditions,  and  if  the  resistance  R 
allows  it,  the  current  will  increase  indefinitely  (i.e.,  limited  by 
the  power  supply)  and  further  transitions  to  a  glow  and  subse¬ 
quently  to  an  arc  discharge  will  be  observed  [6] .  The  duration 
of  the  intermediate  glow  stage  depends  on  a  number  of  param¬ 
eters  (pressure,  electrode  material,  etc.). 

In  the  context  of  this  work,  breakdown  is  defined  as  the 
transition  from  a  nonself-3UStJ>ining  discharge  to  a  Townsend 
discharge  [7] .  The  voltage  across  the  gap  for  this  transition  to 
occur  is  called  the  self -breakdown  voltage  Using  a  similar 

experimental  arrangement,  pulse  experiments  also  have  been 
carried  out.  In  these  experiments,  the  magnitude  of  the  ap¬ 
plied  voltage  pulse  is  less  than  The  purpose  here  is  to 
study  the  space-time  evolution  of  an  initial  bunch  of  electrons 
released  at  the  cathode  by  a  flash  of  ultraviolet  radiation. 
These  experiments  are  referred  to  as  swarm  experiments  [2] . 
Information  obtained  from  both  the  dc  and  the  pulsed  experi¬ 
ments  are  then  used  in  support  of  a  broader  space-time  depen¬ 
dent  theory  (see  Section  III)  and  as  a  tool  to  measure  some  of 
the  fundamental  parameters  associated  with  this  theory. 

In  the  second  approach,  a  coaxial  arrangement  as  shown  in 
Fig.  3  is  used  [8] -[10] .  A  voltage  pulse  of  magnitude  greater 
than  the  static  breakdown  voltage  is  launched  in  the 
coaxial  line  and  it  appears  across  the  electrodes  formed  by  an 
interruption  in  the  center  conductor  of  the  coaxial  line.  The 
voltage  across  the  gap  or  the  current  in  the  transmission  line  is 
recorded  as  a  function  of  time.  A  schematic  indication  of  the 
evolution  of  the  gap  voltage  as  a  function  of  time  is  shown  in 
Fig.  4(a)  and  (b).  Tne  time  elapsed  from  the  application  of 
the  high  voltage  across  the  electrodes  to  the  observed  time  for 
the  collapse  of  the  voltage,  across  the  gap,  i.e,,  To  is  called  the 
observational  time  lag  and  constitutes  what  in  practice  is  called 
the  switch  delay.  A  second  transition  to  a  lower  voltage  may 
be  observed  (i.e.,  in  Fig.  4(a))-  This  corresponds  to  the 
onset  of  an  arc  discharge.  As  the  pressure  of  the  gas  in  the 
gap  and  the  applied  voltage  pulse  is  increased,  the  interval 
-  To  becomes  shorter  in  duration.  Finally,  a  situation 
occurs  where  only  one  transition  is  observed,  as  shown  in 
Fig.  4(b). 

With  this  technique,  it  is  difficult  to  correlate  the  current- 
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Fig.  1 .  Diagram  of  circuit  used  for  the  measurement  of  steady  state 
ionization  currents.  The  gap  electrodes  are  shaped  to  produce  a 
uniform  field  in  the  interelectrode  region.  The  electrode  separation 
is  made  variable  (see  [4]). 


Fig.  2.  Log  IIIq  versus  gap  separation  curves  taken  from  J.  Dutton, 
et  al.  [29]  showing  the  exponential  growth  of  the  field  intensified 
current. 

voltage  measurements  made  with  the  physical  processes  in  the 
gas.  In  particular,  one  cannot  identify  for  certain  the  transi¬ 
tion  occurring  at  Tq.  At  low  pressures  (i.e.,  *^100  torr)  this 
transition  can  be  attributed  to  the  breakdown  of  the  gas,  i.e., 
the  establishment  of  a  self-sustaining  Townsend  discharge 
which  almost  immediately  goes  over  into  a  glow  discharge 
with  maintaining  voltage  (see  Fig.  4(a)).  For  these  condi¬ 
tions,  the  observational  time  lag  consists  of  two  portions;  the 
statistical  time  lag  (the  time  elapsed  from  the  application  of 
the  voltage  to  the  appearance  of  a  breakdown  initiative  elec¬ 
tron  in  the  gap),  and  the  formative  time  lag  (the  time  it  takes 
the  space  charge  to  grow  in  the  gap  which  leads  to  breakdown). 
In  externally  triggered  gas  switches  (laser  or  electron-beam 
triggered,  tor  example)  the  statistical  component  of  the  switch 
delay  is  practically  nonexistent  [11],  [12].  This  type  of 
experiment  represents  a  different  regime  since  the  space  charge 
introduced  by  the  triggering  mechanism  can  be  very  large. 

At  high  pressures  (i.e.,  ^1  atm)  and  for  voltages  a  few  per¬ 
cent  above  (see  Fig.  4(b)),  the  transition  at  Tq  cannot,  with 
certainty,  be  correlated  to  the  breakdown  of  the  gas  as  we 
have  defined  it  [13] .  This  problem  has  been  the  topic  of  con¬ 
siderable  debate  and  will  be  discussed  later  in  Section  IV. 

III.  The  Townsend  Theory  of  Breakdown 

Early  work  in  this  field  was  carried  out  by  Townsend  using 
an  experimental  arrangement  similar  to  that  shown  in  Fig.  I 
[14] ,  In  this  work,  the  current  in  the  external  circuit  as  a 
function  of  electrode  separation  was  found  to  obey  the  follow¬ 
ing  relationship: 

/  =  /oexp(ac/)  (1) 

where  a  is  the  Townsend  primary  ionization  coefficient.  This 
coefficient  corresponds  to  the  number  of  secondary  electrons 
generated  by  a  primar>'  electron  in  moving  1  cm  along  the  field 
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Fig.  3.  Schematic  diagram  of  experimental  arrangement  used  in  pulsed  breakdown  experiments. 


Fig.  4.  Schematic  indication  of  the  evolution  of  gap  voltage  with  time, 
(a)  Low  pressure,  (b)  High  pressure. 

direction.  A  plot  of  In  (///o)  versus  ^  is  a  straight  line  whose 
slope  is  a  as  shown  by  the  lower  portion  of  the  curves  in  Fig.  2, 
As  the  separation  d  increases,  the  curves  begin  to  deviate  from 
this  straight  line  indicating  that  other  processes  begin  to  have 
an  influence  on  the  value  of  the  current.  To  account  for  this 
deviation,  Townsend  proposed  that  the  ions  formed  by  the 
advancing  primary  electrons  contribute  to  the  growth  of  the 
current  by  creating,  as  they  move  towards  the  cathode,  more 
electron-ion  pairs  via  collisions  with  the  neutrals  [14],  A 
number  of  other  “secondary”  processes  have  since  been  found 
which  play  a  larger  role  in  accounting  for  the  faster  than  expo¬ 
nential  growth  of  the  current.  A  number  of  these  processes 
have  been  discussed  by  Llewlyn  Jones  [6] .  The  most  impor¬ 
tant  of  these  secondary  processes  are  ion  impact  on  the  cath¬ 
ode  and  the  photoelectric  effect  at  the  cathode  resulting  from 
the  incident  flux  of  photons  generated  in  the  gap  volume. 


Taking  these  processes  into  account,  the  current  in  the  cir¬ 
cuit  is  found  to  be  given  by  [6] 

U  -  ^  (exp  (ad)  -  l)j 

where  tj  is  a  generalized  Townsend  secondary  ionization  coef¬ 
ficient,  which  depends  on  the  various  secondary  processes 
present  ‘  [6] .  As  shown  in  Fig.  2,  this  expression  describes  the 
experimental  data  quite  well.  As  the  electrode  separation  is 
increased,  note  that  the  slope  tends  toward  infinite.  At  d  = 
the  gap  breaks  down;  i.e,,  /  is  finite  for  Iq  =  0.  The  breakdown 
criterion  is  obtained  from  (2)  by  setting  the  denominator  to 
zero.  This  criterion  is  thus 

^(exp  (oc<i)  -  1)  =  1.  (3) 

a 

As  pointed  out  by  Llewlyn  Jones  [3] ,  this  criterion  for 
breakdown  should  be  taken  more  as  a  physical  interpretation 
of  the  conditions  in  the  gap  than  as  a  derived  mathematical 
criterion,  since  the  theory  as  presented  above  is  a  steady  state 
theory  which  has  no  provisions  for  transient  events  such  as 
breakdown.  Thus  (3)  implies  that  breakdown  occurs  when  the 
secondary  processes  resulting  from  the  motion  of  a  primary 
electron,  and  the  exp  (pcd)  electrons  it  subsequently  creates  in 
traversing  the  gap  (i.e.,  an  avalanche),  are  sufficient  to  regener¬ 
ate  the  primary  electron  at  the  cathode.  In  this  context,  the 
current  becomes  self-sustaining.  Under  transient  conditions, 
the  time  it  takes  for  a  gap  to  breakdown  depends  on  the 
nature  of  the  secondary  process.  Assuming  that  ion  impact  on 
the  cathode  is  the  most  important  mechanism,  the  breakdown 
time  becomes  inversely  proportional  to  the  ion  drift  speed,  i.e., 
in  the  microsecond  range  [6] .  Since  the  processes  are  statisti¬ 
cal,  the  cumulative  effect  of  a  number  of  avalanches  are  needed 
in  order  to  insure  that  the  rejuvenation  process  of  primary 
electrons  will  not  cease,  A  number  of  criteria,  similar  to  (3), 
for  the  breakdown  of  a  gap  have  been  obtained  [7],  [15]. 
They  differ  in  the  functional  form  of  the  coefficient  u,  which 
varies  in  accordance  to  what  is  considered  to  be  the  most 
important  secondary  mechanism. 

Since  ot/p  is  a  unique  function  of  Ejp,  and  further,  if  one 
assumes  that  cu/p  is  also  a  function  of  Ejp  (this  has  been 
experimentally  confirmed  [3]),  then  (3)  represents  a  func¬ 
tional  relationship  between  the  electrode  separation  and  the 
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voltage  at  which  breakdown  occurs.  This  fact  was  first  ob¬ 
served  experimentaUy  by  De  la  Rue  and  MuUer  [16] ,  and  was 
later  extensively  studied  by  Paschen  [17],  Letting  a/p  = 
fi^lp)  2nd  oj/p  =  F(Elp),  (3)  can  be  rewritten  as  [3] 

(exp  (/(£/p)  pd,)  -  1)  =  1 .  (4) 

Since  E  -  V^ld^,  the  above  equation  can,  after  solving  for 
be  rewritten  as 

Vay  =  Hpd,)  (4a) 

which  means  that  the  breakdown  voltage  ^,5  of  the  uniform 
field  gap  is  a  unique  function  <i>  of  the  product  of  pressure  and 
electrode  separation  for  a  given  gas  and  electrode  material. 
This  relationship  is  called  the  Paschen  Law. 

In  certain  regimes,  the  functional  relationship  between  o/p 
andf/p  can  be  approximated  by  [18] 

a/p  =  ^  exp  (-£?/£)  (5) 


nd  (jl0*'®ei»'2; 

Fig.  5.  Comparison  of  the  experimental  and  calculated  values  of  the 
static  breakdown  voltage  as  a  function  of  gap  separation.  The  curve 
is  taken  from  [2,  p.  228] . 


where  A  and  B  are  constants  for  a  given  gas.  Assuming  that 
ion  impact  on  the  cathode  is  the  most  important  secondary 
process,  where  7  is  the  probability  that  a  secondary 

electron  will  be  ejected  from  the  cathode  upon  impact  by  an 
ion  [3] .  Using  (5)  and  the  experimentally  observed  fact  that 
7  is  a  slowly  varying  function  of  £/p  over  a  wide  range,  and 
after  solving  for  ,  (4)  becomes 

V 

*^'C+In(p<fj) 

where 

A 

~l - r-  (6a) 

7^  / 

From  this  equation  it  can  be  seen  (Fig.  5)  that  at  large  values 
of  pdj  the  breakdown  voltage  Fjj,  increases  with  increasing 
pdj.  Similarly,  at  low  values  of  pd^,  increases  with  de¬ 
creasing  pd^.  In  the  transition  from  low  to  large  values  of  pd^, 
the  breakdown  voltage  goes  through  a  minimum  (the  Paschen 
minimum).  This  minimum  breakdown  voltage  is  obtained 
from  (6)  by  differentiating  with  respect  to  pd^  and  equating 
the  result  to  zero.  This  gives  [19] 


pdj  S  200  mmHg  •  cm.  However,  for  large  values  of  pd^,  no 
such  agreement  was  found  in  these  early  experiments. 

At  the  same  time,  experiments  of  Type  II  were  being  con¬ 
ducted  for  large  values  of  pd  and  overvoltage  conditions  [22] . 
Formative  times  of  the  order  of  10''^  s  were  observed.  These 
values  were  too  small  to  be  associated  with  the  7  process  which 
was  considered  to  be  the  most  important  secondary  mecha¬ 
nism  at  the  same  time.  These  results  and  the  deviation  from 
Paschen’s  Law  observed  at  high  pressures,  gave  rise  to  a  num¬ 
ber  of  questions  regarding  the  validity  of  the  Townsend 
Theory  at  large  values  of  pd  [23] ,  [24] .  Three  issues  were  of 
primary  concern.  First,  the  formative  time  lags  in  Type  II 
experiments  were  much  shorter  than  the  microsecond  time  lag 
predicted  from  arguments  compatible  with  the  Townsend 
Theory.  Second,  experiments  of  Type  I,  for  large  values  of  pd, 
seemed  to  indicate  that  the  breakdown  potential  F^j,  was,  to 
a  very  high  degree,  independent  of  the  cathode  material. 
Third,  the  breakdown  was  filamentary  in  character,  which  was 
not  consistent  with  the  Townsend  Theory  of  breakdown.  A 
listing  of  the  objections  to  the  Townsend  Theory  as  result 
of  the  experimental  observations  available  at  the  time  was 
given  by  Meek  [24].  Additional  details  on  the  Townsend 
Theory  can  be  found  in  [3] ,  [6] ,  [7] . 


With  a  voltage  lower  than  ,  it  is  impossible  to  cause  the 

breakdown  of  a  uniform  field  gap  no  matter  what  the  spacing 
or  pressure  is. 

A  number  of  experimental  investigations  (of  Type  I)  to  con¬ 
firm  Paschen’s  Law  have  been  carried  out  [20],  [21].  The 
experimental  results  agreed  quite  well  with  (6)  for  values  of 


IV.  The  Stre-amer  Theory 
As  a  consequence  of  the  above  objections  to  the  Townsend 
Theory,  a  new  theory  was  proposed,  independently,  by  Meek 
[24] ,  Raether  [22] ,  and  Loeb  [25] ,  Meek’s  objectives  were 
to  obtain  a  modified  form  of  (6)  which  would  be  compatible 
with  the  experimental  observations  at  high  values  of  pd.  The 
others  directed  their  attention  to  the  observed  short  time  lags 
and  to  the  filamentary  character  of  the  breakdown. 

In  reality,  theirs  was  more  than  one  theory  [26] .  However, 
the  fundamental  principle  used  by  these  investigators  in  the 
development  of  their  theory  was  common  to  ah.  This  idea  was 
that  at  a  certain  stage  in  the  development  of  a  single  avalanche, 
photoionization  of  the  gas  in  the  interelectrode  space  became 
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the  most  important  mechanism  in  determining  the  breakdown 
of  the  gap.  The  main  difference  between  the  theories  lay  in 
the  stage  at  which  departure  from  the  avalanche-like  develop¬ 
ment  occurred.  For  Meek  and  Loeb,  this  occurred  when  a 
single  avalanche  reached  the  anode,  where  as  for  Raether,  it 
occurred  when  the  avalanche  was  somewhere  in  the  middle  of 
the  gap.  From  this  stage  on,  the  theories  invoked  different 
physical  arguments  to  explain  the  breakdown  of  the  gap. 
However,  they  all  used  photoionization  of  the  gas  surrounding 
the  avalanche  as  the  principal  mechanism. 

According  to  Raether  [23] ,  on  its  way  to  the  anode,  the 
avalanche  reaches  a  critical  dimension  (determined  by  the 
number  of  electrons  in  the  avalanche  or  conversely  by  the 
space  charge  field  produced  by  the  avalanche)  such  that  sec¬ 
ondary  electrons  begin  to  be  generated  just  ahead  of  the  ava¬ 
lanche  by  photoionization  of  the  gas  caused  by  the  ionizing 
radiation  generated  in  the  avalanche.  Of  these  electrons,  there 
is  a  group  located  just  ahead  of  the  avalanche  and  close  to  the 
avalanche  axis.  This  group  is  in  a  region  of  high  field  due  to 
the  enhancement  caused  by  the  avalanche  space  charge.  In 
this  region  of  high  a,  the  electrons  can  multiply  efficiently  (1) 
generating  a  space  charge  cloud  which  rapidly  grows  to  the 
dimension  of  the  parent  cloud,  but  at  a  position  closer  to  the 
anode.  The  process  is  repeated  continuously  to  the  anode 
boundary.  This  progression,  due  to  photoionization,  is  called 
a  streamer  (see  Fig.  6).  Since  the  photons  generating  the  sec¬ 
ondary  electrons  propagate  at  the  speed  of  light,  and  since  the 
subsequent  growth  of  the  space  charge  is  very  fast  (it  occurs  in 
an  enhanced  field),  the  propagation  time  from  the  parent  ava¬ 
lanche  to  the  anode  is  very  small.  Once  the  anode  is  reached, 
Raether  argues,  a  similar  process  begins  to  occur  at  the  cath¬ 
ode  end  of  the  parent  avalanche.  There  the  photoelectrons 
generated  are  accelerated  towards  the  avalanche  extending  the 
ion  sheath  of  the  parent  avalanche  towards  the  cathode 
(Fig.  6).  Breakdown  occurs  immediately  upon  the  space 
charge  cloud  reaching  the  cathode. 

On  the  other  hand,  Loeb  and  Meek  considered  only  the 
cathode  directed  streamer.  That  is,  they  assumed  that  the 
parent  avalanche  becomes  “critical”  just  in  front  of  the  anode. 

What  determines  the  critical  stage?  Meek  arbitrarily  chose 
it  to  be  when  the  field  due  to  the  avalanche  is  equal  to  the 
external  field.  Using  this  argument  he  found  that  the  break¬ 
down  equation  (i.e.,  Paschen’s  Law)  becomes  [24] 

oc  El 

(alp)  (pd)  +  In  —  =  14.46  +  In  —  -  In  pd  +  In  (9) 

P  P  2 

Raether,  on  the  other  hand,  defined  the  critical  stage  from 
experimental  observations  of  avalanche  development  in  a 
cloud  chamber.  He  observed  that  streamers  developed  when 
[23] 

aX^,c:i20 

where  is  the  average  position  of  the  avalanche  when  it 
becomes  critical.  Using  this,  he  calculated  a  formative  time 
for  breakdown.  He  argued  that  since  the  streamer  propagation 
velocity  is  high  (of  the  order  of  the  speed  of  light),  the  forma- 
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Fig.  6.  Sketch  of  the  propagation  of  a  streamer  due  to  gas  ionizing 
radiation,  (a)  Anode  directed  streamer,  (b)  Cathode  directed 
streamer. 


tive  time  is  basically  the  time  it  takes  an  avalanche  to  become 
critical  fcr-  This  time  is  given  by  [23] 


(10) 


where  is  the  electron  drift  velocity  in  the  applied  field. 
The  values  obtained  from  (9)  and  (10)  were  very  close  to  ex¬ 
perimental  observations. 

Further  support  for  the  streamer  theory  came  from  Type  II 
experiments  by  Fletcher  [8].  He  also  tried  to  improve  the 
mathematical  framework  of  the  theory;  however,  the  funda¬ 
mental  premise,  i.e.,  photoionization  of  the  gas  was  left  un¬ 
questioned.  Similar  types  of  experiments  were  carried  out  by 
Felsenthal  and  Proud  [28] .  They  were  able  to  explain  their 
results,  however,  using  a  modified  pulsed-microwave  break¬ 
down  theory  [28].  Further  details  on  the  streamer  theory 
may  be  found  in  [8] ,  [23] -[25] ,  [28] . 


V.  The  Poststreamer  Era 

In  the  preceding  two  sections,  the  two  principal  models 
proposed  to  explain  the  initial  phase  of  breakdown  were  dis¬ 
cussed.  At  low  values  of  pd  (i,e.,  pd<  200  mmHg  *  cm),  it 
was  generally  agreed  that  the  Townsend  Theory  gave  a  satis¬ 
factory  explanation  of  the  observed  phenomena.  At  the  high 
values  of  pd,  there  was  no  general  agreement  as  to  a  theory  for 
breakdown  [3] ,  [7] . 

At  this  time,  two  schools  regarding  breakdown  at  high  values 
of  pd  developed.  One  group  contended  that  the  Townsend 
mechanism,  which  involves  the  interaction  of  primary  and 
secondary  processes  associated  with  the  cathode,  was  still 
operative.  The  other  group  supported  the  streamer  mecha¬ 
nism,  The  burden  of  proof  lay  with  the  proponents  of  the 
Townsend  mechanism  since  the  inconsistencies  between  the 
Townsend  theory  and  experiments  were  still  unresolved.  The 
work  by  this  later  group  concentrated  in  three  areas. 

Firstly,  experiments  of  Type  I  were  carried  out  with  increas¬ 
ingly  better  accuracy  by  Llewelyn  Jones  and  his  colleagues  at 
Swansea  [3]-[5].  The  purpose  of  these  experiments  was  to 
show  that  the  apparent  departure  from  Paschen’s  Law  ob¬ 
served  at  high  pressures  was  due  to  poor  experimental  condi¬ 
tions.  Llewelyn  Jones  argued  that,  at  high  pressures,  a  very 
stable  voltage  source  was  necessary  in  order  to  observe  the 
effects  of  cathode  material  on  the  breakdown  voltage  (i.e.,  the 
up  curving  in  Fig.  2).  This  may  be  seen  from  (4).  The  frac- 
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tional  change  in  breakdown  voltage  is  given  by  [3] 
~  [pdj: 


(11) 


For  large  values  of  is  approximately  constant, 

so  that  AKjb/Fsb  ^  It  follows  that  for  large 

can  become  small  enough  to  be  unobserved  with  the 
experimental  voltage  supplies  used  in  earlier  experiments.  In 
the  last  measurements  done  at  Swansea,  the  effects  of  second- 
ary  processes  were  observed  for  pd^  values  up  to  12  000 
mmHg  •  cm  [29] . 

Secondly,  a  space/ time  dependent  formulation  of  break¬ 
down  based  on  the  Townsend-type  mechanism  taking  the 
space  charge  and  the  photoelectric  effect  into  consideration 
was  developed.  The  idea  was  to  show  that  short  breakdown 
times  (^10*^  s)  could  be  obtained  under  a  modified  Town¬ 
send  formulation.  The  first  attempt  in  this  direction  was  by 
Davidson  [30]  who  obtained  an  analytic  solution  for  the 
growth  of  ionization  in  the  gap  assuming  that  space  charge 
effects  were  negligible.  The  inclusion  of  space  charge  effects 
made  the  solution  to  the  problem  analytically  untractable:  a 
numerical  solution  of  the  system  of  equations  [31]  was  re¬ 
quired.  This  study  showed  that  short  time  lags  ('^10*'^  s) 
could  indeed  be  obtained  using  a  Townsend  model.  Moreover, 
it  was  shown  by  Dickey  [13]  that  even  the  very  short  time 
lags  observed  by  Fletcher  [8]  (^10"^  s)  could  be  explained 
if  the  effect  of  the  motion  of  the  space  charge  in  the  gap  on 
the  external  circuit  is  taken  into  consideration.  He  noted  that 
if  a  large  space  charge  is  moving  across  the  gap,  a  large  current 
can  be  generated  in  the  external  circuit  which  can  (because  of 
the  drop  in  voltage  across  the  external  resistance)  cause  the 
collapse  of  the  voltage  across  the  gap.  Thus  a  voltage  collapse 
is  observed  even  before  the  gap  is  bridged  by  the  space  charge. 
The  values  he  obtained  for  the  time  lags  using  this  argument 
were  close  to  those  observed  by  Fletcher  [13].  Thus  what 
constitutes  “breakdown”  under  these  conditions  is  not  only 
different  from  the  definition  given  earlier,  but  is  difficult  to 
define. 

The  third  area  was  the  issue  of  interpreting  luminosity  obser¬ 
vations.  Raether  first  noted  from  luminosity  observations  that 
the  velocity  of  propagation  of  a  streamer  was  of  the  order  of 
10®  cm/s  [23].  Subsequently,  a  number  of  photographic 
experiments  were  carried  out  which  yielded  conflicting  results 
as  far  as  the  development  of  the  streamer  [32]-[34] .  In  ex¬ 
periments  by  Wagner  [32] ,  the  streamer  was  observed  to 
have  two  different  velocity  stages;  whereas,  experiments  by 
Chalmers*  [33]  showed  there  is  a  single  propagating  velocity. 
These  discrepancies  may  be  due  to  sensitivity  differences  in 
the  streak  camera  or  film  [34] .  However,  as  was  pointed  out 
by  Llewelyn  Jones  [6],  luminosity  measurements  could  be 
very  misleading.  He  disagreed  with  Raether’s  conclusions 
about  streamer  propagation  obtained  from  the  luminosity 
observations.  He  suggested  that  the  “streamer”  seen  by 
Raether  was  merely  a  luminous  “phase”  effect  and  that  no 
real  physical  process  was  developing  in  the  gap  [6] .  That  is, 
that  the  propagation  of  the  luminosity  front  was  due  to  the 


temporal  growth  of  the  space  charge  to  a  level  which  can  be 
recorded  from  its  luminosity,  and  not  to  the  spatial  propaga¬ 
tion  of  the  space  charge.  Computer  simulations  subsequently 
were  carried  out  at  Swansea  to  show  that  this  interpretation 
was  indeed  correct  [35] ,  [36] .  The  results  seem  to  agree 
quite  well  with  an  experiment  by  Wagner  [32] .  The  simula¬ 
tions  did  not  take  photoionization  into  account.  A  compre¬ 
hensive  compilation  of  this  work  is  given  in  [37] .  Reference 
[38]  gives  a  review  of  luminosity  studies. 

Since  the  time  lags  predicted  from  the  Townsend  and 
streamer  theories  were  different,  a  discontinuity  in  the  time 
lag  versus  pd  curve  should  be  observed,  indicating  the  transi¬ 
tion  from  one  regime  to  the  other.  The  first  experiment  (Type 
II)  designed  to  look  for  such  transition  was  carried  out  by 
Fischer  and  Bederson  [39] .  The  results  were  rather  dramatic 
and  completedly  changed  the  traditional  boundary  separating 
the  two  mechanisms.  They  found  that  for  applied  voltages 
slightly  above  long  time  lags  were  observed  forpcf  values 
up  to  1000  mmHg  •  cm.  The  time  lags  decreased  continuously 
as  a  function  of  the  percent  overvoltage.  This  decrease  they 
attributed  to  the  role  played  by  the  avalanche  space  charge 
and  by  the  photo-electric  effect.  It  was  concluded  from  these 
experiments  that  percent  overvoltage,  and  not  pd,  was  the 
principal  parameter  determining  the  value  of  time  lag.  This 
was  later  confirmed  by  the  computer  analysis  of  Ward  [40] . 
Allen  and  Phillip  [41]  later  extended  these  experiments  to 
larger  values  of  percent  overvoltage  and  to  a  number  of  differ¬ 
ent  gases.  Their  results  showed  a  change  in  the  slope  of  the 
time  lag  versus  percent  overvoltage  curve  for  some  of  the  gases 
studied.  This  they  interpreted  as  an  indication  of  a  transition 
from  the  Townsend  regime  to  the  streamer  regime. 

It  should  be  mentioned  that  a  large  amount  of  work  in 
atomic  physics  was  carried  out  during  this  period  which  helped 
elucidate  the  importance  of  the  various  processes  occurring 
during  breakdown.  Moreover,  they  provided  increasingly  more 
accurate  values  for  the  various  transport  coefficients  used  in 
the  numerical  calculations.  For  more  information  on  this  type 
of  experiments  the  reader  should  consult  [2] . 

For  the  streamer  group,  a  number  of  experiments  were 
carried  out  to  measure  the  spectrum  of  the  radiation  generated 
by  the  avalanche  and  the  absorption  coefficient  for  the  “ioniz¬ 
ing”  radiation  [43] ,  [44] .  During  this  time,  very  little  was 
added  to  the  fundamental  theories  of  Raether,  Meek,  and 
Loeb,  until  Lozanskii  [15]  in  1968  tried  to  explain  the  nature 
of  the  process  for  photoionization  ahead  of  the  avalanche. 
This  has  been  an  ill-defined  issue  in  the  streamer  theory  and 
has  been  questioned  by  many  authors  [3],  [45].  The  main 
issue  has  been  whether  or  not  photons  of  high  enough  energy 
to  photoionize  are  produced  by  the  avalanche,  and  if  they  are 
produced,  can  they  escape  the  avalanche?  This  question  is 
particularly  difficult  to  answer  in  experiments  with  relatively 
pure  gases  [27] .  Lozanskii  proposed  that  associative  ioniza¬ 
tion  of  atoms  excited  by  photons  of  energy  lower  than  that 
required  to  ionize  are  responsible  for  the  formation  and  propa¬ 
gation  of  the  streamer.  Assuming  a  Lorentzian  profile  for  the 
spectral  line  of  interest,  he  obtained  an  expression  for  the 
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absorption  coefficient  of  these  lower  energy  photons  [15], 
Because  of  the  spectral  line  width,  the  off  resonant  photons 
can  escape  the  avalanche  and  excite  atoms  just  ahead  of  the 
avalanche.  He  then  estimated  the  number  of  electrons  that 
will  be  produced  via  the  chemical  reaction 

+  (12) 

where  is  an  excited  atom  or  molecule  and  B  is  an  atom  or 
molecule  in  the  ground  state.  The  cross  section  for  this  reac¬ 
tion  is  rather  large.  In  helium,  for  example,  the  cross  section 
for  the  reaction  [27] 

He*  +  He->He^ +e“  (13) 

has  a  value  of  cm^.  The  number  of  electrons  pro¬ 

duced  ahead  of  an  avalanche  via  this  process  is  sufficient  to 
create  and  maintain  a  streamer.  However,  the  occurrence  of 
this  (or  any  other  process)  has  not  been  confirmed  experi¬ 
mentally,  Computer  simulations  for  overvolted  gaps  (<:  20- 
percent  self-breakdown)  have  been  carried  out  by  Kline  [46] 
and  by  Yoshida  and  Tagashira  [47] .  The  absorption  coeffi¬ 
cient  used  by  Yoshida,  et  ai,  for  the  radiation  was  consistent 
with  the  Lozanskii  mechanism.  Good  agreement  was  obtained 
between  these  simulations  and  experiments  done  by  Wagner 
[32]  and  Koppitz  [48]  for  overvolted  gaps.  In  Koppitz’s 
experiments,  the  breakdown  was  large  in  cross  section  and 
thus  it  was  more  justifiable  to  apply  the  one-dimensional 
theory  used  in  the  computations.  The  large  discharge  cross 
section  was  achieved  by  illuminating  the  cathode  surface  with 
ultraviolet  radiation.  One  should  recall  that  success  in  model¬ 
ing  Wagner’s  experiment  was  also  achieved  without  taking 
photoionization  into  account  [35] . 

It  follows  from  the  works  mentioned  above  that  a  properly 
formulated  theory  (i.e.,  including  space  charge  and  photoelec¬ 
tric  effects)  based  on  the  Townsend  mechanism  is  able  to  ade¬ 
quately  describe  the  breakdown  processes  of  a  gap  for  large 
values  of  pd,  as  well  as,  for  applied  voltages  of  a  few  percent 
above  self-breakdown  (up  to  20  percent);  and  that  the  early 
objections  to  the  theory  were  due  to  experimental  “idiosyn- 
crases.”  Thus  the  regime  of  validity  of  the  two  models  was 
apparently  established:  Townsend  theory  for  low  overvoltage 
(<20  percent),  and  streamer  theory  at  high  overvoltage 
(^  20  percent). 

At  very  high  overvoltages  (two  to  three  times  self¬ 
breakdown),  the  character  of  the  breakdown  once  again 
changes  [49] :  a  broad  channel  is  observed  as  opposed  to  a 
filamentary  channel.  This  has  been  explained  by  Babich  and 
Stankevich  [50]  in  terms  of  runaway  electrons  which  can  be 
produced  in  the  gap  at  these  high  fields.  Experimental  evi¬ 
dence  of  these  runaway  electrons  exists;  however,  these  results 
are  not  time  resolved  so  that  it  is  not  known  when  during  the 
pulsed  discharge  the  runaway  electrons  were  produced  [51] . 

For  voltages  below  those  of  the  above  experiments,  Mesyats 
and  colleagues  [52]-[54]  obtained  values  for  the  formative 
time  which  were  not  in  agreement  with  those  obtained  using 
the  streamer  model.  They  noted  that  the  “breakdown  time” 
depended  on  the  number  of  initiating  electrons  at  the  cath¬ 
ode.  For  large  numbers  of  initiatory  electrons  ("^10"^),  Mesyats 
derived  an  expression  for  the  “breakdown  time”  using  argu¬ 


Fig.  7.  Schematic  indication  of  axial  electric  field  in  avalanche  region 

[55]. 

ments  similar  to  Dickey’s  [53] .  He  pointed  out  that  this  was 
indeed  the  initial  condition  in  Fletcher’s  experiment.  For  low 
number  of  initiatory  electrons  (which  he  called  single  electron 
initiation),  the  “breakdown  time”  was  long  compared  to  that 
obtained  from  the  streamer  model  [54] .  He  argued  that  since 
the  time  for  an  avalanche  to  become  critical,  at  these  high 
values  of  electric  field,  was  much  shorter  than  the  average  life¬ 
time  of  the  excited  atoms,  photoionization  played  no  role  in 
the  breakdown  process;  moreover,  that  the  avalanche  stops 
growing  due  to  the  decrease  in  field  at  the  avalanche  head  due 
to  the  local  space  charge.  At  this  point,  a  new  avalanche 
develops  ahead  of  the  parent  avalanche  which  goes  through 
the  same  cycle.  Thus  an  avalanche  chain  is  formed  which 
slowly  drifts  across  the  gap.  He  called  this  the  avalanche  chain 
model,  where  a  linear  chain  of  avalanches,  originating  at  the 
cathode  and  supported  by  the  photoelectric  effect  at  the  cath¬ 
ode,  is  responsible  for  bridging  the  gap  [54] . 

The  present  state  of  affairs  is  that  the  domain  of  overvoltage 
is  divided  as  follows:  Townsend  model,  streamer  model,  ava¬ 
lanche  chain  model,  and  the  continuous  acceleration  model. 

Recently,  Kunhardt  and  Byszewski  [55]  introduced  a  model 
for  the  development  of  breakdown  above  the  Townsend 
regime.  The  model  offers  the  attractive  feature  of  the  unifica¬ 
tion  of  all  the  proposed  breakdown  phenomena.  That  is,  it 
gives  a  continuous  picture  of  the  breakdown  above  the  Town¬ 
send  regime,  and  merges  into  the  Townsend  avalanche  picture 
as  the  voltage  is  reduced. 

In  this  model,  the  energy  distribution  function  for  electrons 
in  the  advancing  avalanche  is  assumed  to  have  two  compo¬ 
nents:  fast  electrons  and  slow  (thermal)  electrons.  The  fast 
electrons  can  “run  away”  from  the  avalanche.  This  happens 
because  the  effective  retarding  force  on  an  electron  moving 
through  a  neutral  gas  decreases  with  increasing  velocity,  in  the 
case  of  electrons  possessing  sufficiently  high  energy  (i.e., 
u  ^  3-5 €/,  where  6/  is  the  ionization  energy).  The  energy 
threshold  for  these  runaway  electrons  is  determined  by  the 
electric  field  strength  (Fig.  7(a))  and  the  gap  parameters.  Once 
these  fast  electrons  leave  the  avalanche,  most  of  them  no 
longer  meet  the  runaway  condition  and  become  “trapped.” 
This  is  due  to  the  fact  that  they  enter  a  region  of  decreasing 
field  ahead  of  the  avalanche,  and  the  energy  they  gain  from 
the  field  along  their  trajectory  is  not  enough  to  overcome  the 
losses  (Fig.  8).  The  “trapping”  distance  of  these  electrons  is  a 
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Fig.  8.  Retarding  force  versus  electron  energy  [55] . 
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Fig.  9.  Election  energy  versus  distance  from  avalanche  front  [55  ] . 


function  of  their  initial  energy  (Fig.  9).  Thus  these  “seed” 
electrons,  which  are  continuously  “emitted”  from  the  ava¬ 
lanche,  multiply  at  various  distances  from  the  parent  avalanche 
rapidly  extending  the  avalanche  space  charge  towards  the 
anode.  On  the  cathode  side,  the  photoelectric  effect  is  as¬ 
sumed  to  be  the  primary  mechanism  for  generating  secondary 
electrons  which  are  subsequently  accelerated  towards  the  high 
field  region  of  the  parent  avalanche. 

This  model  has  been  used  to  explain,  among  other  things, 
the  variations  in  the  time  lag  versus  percent  overvoltage 
observed  by  Men  and  Phillips  [41]  for  different  gases,  and 
the  “broad”  breakdown  observed  by  Koppitz  [48]  when  the 
initial  ultraviolet  irradiation  covered  a  large  portion  of  the 
cathode  surface. 

The  most  important  concept  associated  with  an  avalanche  is 
that  of  the  nature  of  the  electron  distribution  function  [56] . 


At  present,  this  is  one  of  the  principal  areas  of  investigation 
[57]-[60].  Because  of  the  complexity  of  the  problem,  the 
work  is  mainly  computational,  although  some  analytical  at¬ 
tempts  have  been  made  [61] .  Computationally,  Monte  Carlo 
[57]  ,  [58]  and  Boltzman  [59],  [60]  codes  with  varying 
degrees  of  sophistication  have  been  used.  The  thrust  is  to 
arrive  at  the  distribution  function  for  electrons  in  regions  of 
high  electric  fields  and  large  spatial  gradients.  No  model  can 
conclusively  be  validated  until  this  distribution  is  obtained. 

Our  ability  to  manipulate  the  switch  delay  time  will  depend 
on  our  understanding  of  the  evolution  of  the  electron  energy 
distribution  function  under  these  highly  “nonequilibrium” 
conditions. 
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A  model  for  the  development  of  electrical  breakdown  in  dense  gases  is  presented.  It  describes  the  imtial 
phase  of  breakdown  in  the  regime  where  the  Townsend  avalanche  mechanism  does  not  apply.  The  main 
features  of  the  model  are  as  follows:  (1)  It  gives  a  continuous  picture  of  the  development  both  in  the 
structure  of  the  breakdown  and  the  physics  of  the  processes,  and  (2)  it  is  based  on  electron  kinetics,  so  that 
the  theory  is  general  in  scope.  In  Ught  of  this  model  a  brief  discussion  of  experimental  results  is  given. 


I.  INTRODUCTION 

Two  basic  models  have  been  proposed  to  explain 
the  initial  phase  of  electrical  breakdown  in  gases 
at  hi^-  pressures.  These  are  commonly  known  as 
the  Townsend  avalanche  model^  and  the  streamer 
model.*  The  mechanism  playing  a  role  in  the  first 
model  involves  the  interaction  of  primary  and 
secondary  processes  of  which  cathode  processes 
are  generally  found  to  predominate  in  uniform 
fields.*  Photoionization  of  the  gas  In  the  inter- 
electrode  volume  forms  the  basis  for  the  second 
model.*  Over  die  last  30  years,  considerable  ef¬ 
fort  has  gone  into  establishing  the  regime  of  valid¬ 
ity  of  both  models.  Experimental  evidence  was 
sought  for  the  transition  between  the  two  regimes, 
such  as  an  abrupt  decrease  in  the  formative  time 
of  the  breakdown  as  the  voltage  was  increased 
above  but  near  the  self-breakdown  voltage.  The 
results  were  not  conclusive.^ 

The  Townsend  avalanche  model  was  initially 
thought  to  apply  only  at  low  pressures  and  for 
voltages  near  self-breakdown.  This  regime,  how¬ 
ever,  has  recently  been  expanded  by  the  experi¬ 
mental  and  theoretical  works  at  Swansea*’*  and  by 
the  computer  simulations  of  Ward,*^  to  include  the 
high-pressure,  small  overvoltage  regime  (^20%). 
This  was  the  regime  for  which  the  streamer  model 
was  originally  proposed.  Unlike  the  Townsend 
avalanche  theory,  the  foundations  of  the  streamer 
theory  have  been  very  vague  from  the  start^  It  is 
important  to  note  that  the  “  streamer  moder  is 
actually  more  than  one  model,  each  based  on  dif¬ 
ferent  hypotheses.*  The  clearest  physical  picture 
of  the  processes  associated  with  a  streamer  model 
has  been  given  by  Lozanskii,®  who  proposed  that 
associative  ionization  of  atoms  excited  by  photons 
of  lower  energy  than  those  required  to  ionize,  are 
responsible  for  the  formation  and  propagation  of 
the  streamer.  The  initial  support  for  the  streamer 
model  came  from  experiments  by  Rogowski^®  and 
Raether,*  and  later  by  Fletcher, “  who  measured 
values  for  the  formation  time  somewhat  in  agree¬ 


ment  with  those  obtained  theoretically.  Lumino¬ 
sity  measurements  initially  used  by  many  workers 
as  further  evidence  in  support  of  the  theory,  have 
been  shown,  in  certain  cases,  to  be  very  mis¬ 
leading.^*  Moreover,  as  pointed  out  by  Mesyats,^* 
agreement  in  the  value  of  the  formative  time  and, 
we  may  add,  in  the  structure  of  the  breakdown 
process,  is  a  necessary  but  not  a  sufficient  crite¬ 
rion  for  the  validity  of  the  theoretical  premises. 
Agreement  in  the  formation  time  was  also  obtained 
by  Dickey,^'*  although  his  theory  was  based  on  a 
completely  different  hypothesis. 

At  very  high  overvoltages  (two  to  three  times 
self -breakdown  voltage),  Stankevich  and  Kalinin^* 
have  observed  that  the  structure  of  the  breakdown, 
once  again,  changes;  i.e.,  a  broad  channel  break¬ 
down  is  observed  as  opposed  to  a  filamentary 
channel.  This  has  been  explained  by  Babich  and 
Stankevich^®  in  terms  of  runaway  electrons  from 
the  main  avalanche.  We  must  point  out  that  there 
are  a  few  inconsistencies  in  the  paper  of  Babich 
et  aL  regarding  the  criterion  for  the  existence  of 
these  runaway  electrons.  Experimental  evidence 
of  these  runaway  electrons  exists.^"'  These  experi¬ 
ments,  however,  were  not  time  resolved  so  that 
it  is  not  known  when  during  the  breakdown  runaway 
electrons  were  produced. 

At  voltages  below  this  regime,  Mesyats^®  obtained 
values  for  the  formative  time  which  were  not  in 
agreement  with  those  obtained  using  the  streamer 
model.  He  then  proposed  the  avalanche  chain 
model,  where  a  chain  of  avalanches  originating 
at  the  cathode  and  supported  by  the  photoelectric 
effect  at  the  cathode  is  responsible  for  bridging 
the  gap. 

Thus,  at  present,  the  streamer  model  is  thought 
to  qualitatively  describe  the  breakdown  processes 
at  high  pressures  and  overvoltages  (above  20% 
self-breakdown),  butbelow  Mesyats’ avalanche  chain 
regime. 

In  this  paper,  we  present  a  model  to  describe  the 
breakdown  processes  lot  all  voltages  and  pressures 
above  the  Townsend  avalanche  domain.  We  call  it  the 
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two- group  model  due  to  the  similarity  between  this 
problem  and  the  neutron  transport  problem  in  nuclear 
reactor  physics*  It  isbased  on  electron  kinetics,  and 
it  is  thus  applicable  to  abroad  range  of  conditions 
(breakdown  in  pure  gases,  for  example).  The  model 
offers  the  attractive  feature  of  the  unification  of  all  the 
proposed  breakdown  phenomena;  that  Is,  it  provides 
a  continuous  picture  of  the  development  of  the  initial 
phase  of  breakdown  above  the  Townsend  regime  both 
in  the  structure  of  the  breakdown  and  the  physics  of 
the  processes.  It  merges  into  the  T  owns  end  avalanche 
picture  as  the  voltage  is  reduced. 

In  Bee.  n,  we  will  present  the  model,  devoting 
our  attention  to  the  physical  processes  that  occur 
and  that  lead  to  the  formation  of  breakdown.  We 
then  proceed,  in  Sec,  HI,  to  attempt  to  formulate 
a  mathematical  framework  for  the  theory.  In  Sec. 
nr,  we  discuss  some  experimental  results  in  light 
of  the  theory  presented. 

n.  TWO-GROUF  MODEL  OF  BRE.4KDOWN 

We  will  develop  a  physical  picture  of  the  elec¬ 
trical  breakdown  of  a  gas  at  high  pressures  and 
for  voltages  so  that  the  effects  of  the  produced 
space  charge  cannot  be  neglected.  This  regime 
lies  above  the  Townsend  avalanche  regime.  It  is 
our  purpose  in  this  section  to  give  a  clear  physical 
picture  of  the  processes  responsible  for  increasing 
die  propagation  speed  of  the  avalanche  towards 
the  anode  and  for  causing  the  filamentary  appear¬ 
ance  of  the  evolving  avalanche. 

Oir  model  is  based  on  electron  kinetics  rather 
than  photon  induced  processes.  Electrons  with  a 
broad  spectrum  of  energies  are  always  present  in 
an  avalanche.  On  the  other  hand,  the  existence  of 
photons  with  enough  energy  and  range  to  photo- 
ionize  the  gas  outside  the  avalanche  has  always 
been  a  questionable  matter,  especially  in  relative¬ 
ly  pure  gzsesJ  The  possibility  of  multiphoton 
ionization  would  require  photon  fluxes  not  available 
in  the  avalanche  stage. 

A  sketch  of  the  experimental  situation  we  want 
to  consider  is  shown  in  Fig.  1.  Also  shown  is  an 
avalanche  propagating  towards  the  anode.  The 
shape  of  the  avalanche  is  indicative  of  the  pro¬ 
cesses  at  play  in  the  development  up  to  its  present 
location  in  the  gap.  The  avalanche  is  divided  into 
three  stages  (see  Fig,  1).  This  division  is  strictly 
for  illustrative  purposes,  since  a  continuous 
transformation  is  actually  taking  place.  In  stage 
I,  the  radial  dimensions  of  the  avalanche  is  deter¬ 
mined  primarily  by  diffusion  processes.  The 
avalanche  radius  is  given  by^ 

'  U) 

where  D  is  the  electron  diffusion  coefficient.  In 
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FIG,  1.  Sketch  of  the  developed  avalanche  in  the 
cathode-anode  volume. 

the  regime  of  interest  (i.e,,  for  voltages  >20'^ 
self-breakdown),  the  time  of  development  is  very 
short  so  that  very  little  expansion  occurs.  We  can 
estimate  the  value  of  the  radius  for  the  case  of 
breakdown  in  nitrogen,  at  atmospheric  pressure. 
Here,  cmVsec,  sec,  and  from  (1), 

r^-7.2x  10'^  cm. 

As  the  number  of  electrons  in  the  avalanche  in¬ 
creases,  electrostatic  repulsion  begins  to  play  a 
role  in  the  expansion  of  the  avalanche.  In  this  re¬ 
gime  (stage  n),  the  avalanche  radius  increases 
exponentially  with  time.  This  radius  may  be  esti¬ 
mated  using  the  expression^® 


where  e  is  the  electron  charge,  is  the  external 
field,  a  is  the  Townsend  primary  ionization  co¬ 
efficient,  z  is  the  average  position  along  the  gap 
axis  of  the  avalanche  with  z  =  0  being  the  cathode 
position,  and  is  the  permittivity  of  free  space. 
For  breakdown  in  nitrogen  at  1  atm  and  an 
applied  field  of  40  kV/cm  (Le.,  -“25%  above  self¬ 
breakdown),  for  which  af'56  cm*S  the  radius  of 
the  avalanche  at  2  =  0,05  cm  is  r,  =  0.14  cm.  The 
radial  expansion  is  somewhat  slowed  down  in  stage 
ni  as  space-charge  neutralization  begins  to  take 
place. 

Throughout  these  three  stages,  as  the  electron 
and  ion  densities  in  the  avalanche  increase,  two 
things  are  happening:  First,  a  highly  nonequili- 
brium  electron  energy  distribution  is  evolving^® 
and  second,  the  external  field  is  becoming  dis¬ 
torted  due  to  the  presence  of  the  space  charge.  It 
is  the  combination  of  these  two  effects  that  play 
the  major  role  in  the  following  stage,  'Ae,,  stage 
rV.  A  qualitative  picture  of  the  longitudinal  field 
at  the  axis  of  the  avalanche  as  a  function  of  dis¬ 
tance  and  just  prior  to  the  time  the  avalanche 
“enters"'  into  stage  IV  is  given  in  Fig.  2.  The  r 
=  0  plane  corresponds  to  the  front  of  the  avalanche. 
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FIG.  2.  Sketch  of  the  magnitude  of  the  longitudinal 
electric  field  along  the  axis  of  the  avalanche. 

Here,  the  field  attains  the  maximum  average  value 
it  can  have  in  the  gap  and  it  decreases  away 
from  this  plane.  The  rate  of  decrease  is  slower 
forx>0. 

The  electron  energy  distribution  in  the  avalanche 
is  far  from  equilibrium.  It  is  enriched  with  high- 
energy  electrons,  a  consequence  of  the  large  elec¬ 
tric  field.  Furthermore,  it  is  anisotropic  in  the 
high-energy  region;  it  acquires  a  directed  charac¬ 
ter  along  the  field  (this  will  be  further  discussed). 
For  the  purpose  of  discussion,  we  may  think  of 
the  electron  energy  distribution  as  being  formed 
of  two  groups  of  electrons:  the  "main”  electron 
distribution  and  the  “fast”  electrons. 

It  is  possible  for  some  of  these  fact  electrons  to 
become  runaways,  that  is,  they  may  continuously 
gain  energy  from  the  field.  This  happens  because 
the  effective  retarding  force  on  an  electron  moving 
through  a  neutral  gas  decreases  with  increasing 
velocity,  in  the  case  of  electrons  possessing  a 
sufficiently  high  energy  (i.e.,  where 

is  the  ionization  energy).  The  energy  threshold 
for  electron  runaway  is  determined  by  the  magni¬ 
tude  of  the  electric  field.  The  larger  the  field, 
the  lower  the  threshold  energy  and  consequently 
the  larger  number  of  electrons  can  run  away. 

In  Fig.  2,  the  region  where  electrons  are  most 
likely  to  run  away  is  labeled  the  ‘'injection  region.” 
In  this  region  of  high  electric  field,  the  runaway 
threshold  energy  is  the  lowest  and  is  determined 
by  the  maximum  field  intensity  Since  the  dis¬ 
tribution  is  anisotropic  in  this  energy  regime,  the 
runaway  electrons  are  accelerated  out  of  this  re¬ 
gion  and  are  injected  into  the  region  ahead  of  the 
avalanche  where  the  field  is  decreasing.  The 
number  of  particles  injected  is  maximum  at  the 
axis.  Here,  the  space-charge  fields  and  the  ex¬ 
ternal  field  are  colinear,  thus  giving  the  maximum 
total  field,  hence,  the  lowest  injection  energy.  As 
we  move  away  from  the  axis,  the  space-charge 
field  and  the  external  field  are  no  longer  collinear, 
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FIG.  3.  Development  of  a  filamentary  channel  ahead 
of  the  avalanche. 

so  that  the  total  field  is  reduced,  thus  increasing 
the  runaway  threshold  energy.  We  thus  observe 
that  there  is  an  injection  cone  with  a  maximum  on 
the  axis  [see  Fig.  3(a)l. 

Once  injected,  most  of  the  fast  electrons  no 
longer  meet  the  runaway  conditions  and  become 
“trapped”;  that  is,  the  energy  they  gain  along 
their  trajectory  is  not  enough  to  overcome  the 
losses  (see  Sec.  III).  The  trapping  distance  (i.e., 
the  distance  from  the  avalanche  front  to  where  they 
become  trapped)  is  for  a  given  gas  a  function  of 
the  initial  injection  energy  of  the  electrons  and  the 
slope  of  the  electric  field  ahead  of  the  avalanche. 

Thus,  just  prior  to  stage  IV,  we  may  think  of  the 
avalanche  as  a  localized  distribution  of  electrons. 
The  beginning  of  stage  IV  [Fig.  3(a)]  is  marked  by 
the  “burst”  of  the  avalanche  along  its  axis,  followed 
by  the  ejection  of  high-velocity,  electrons  and 
their  subsequent  capture  at  varying  distance  from 
the  original  avalanche  position  where  they  start  to 
generate  avalanches  of  their  own. 

The  captured  electrons  ionize  the  gas  extend¬ 
ing  the  boundary  of  the  main  avalanche  along 
a  filamentary  channel  centered  at  the  axis  of  the 
avalanche  [see  Fig.  3(b)].  Once  the  channel  starts 
to  narrow,  the  electrical  field  just  ahead  of  the  tip 
increases.  Fast  electrons  are  continuously  in¬ 
jected  into  the  region  ahead  of  the  tip  and  again 
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are  trapped.  This  process  accelerates  the  devel¬ 
opment  of  the  avalanche  tip  towards  the  anode. 

This  was  shown  by  Raether,’  although  his  ex¬ 
planation  was  based  on  the  electrons  being  pro¬ 
duced  by  photoionization  ahead  of  the  avalanche 
rather  than  by  captured  runaway  electrons.  Thus 
a  filamentary  channel  (note  that  the  channel  is 
maintained  narrow  due  to  the  injection  cone  pheno¬ 
menon)  evolves  from  the  main  avalanche  whose 
structure  is  similar  to  the  classical  streamer, 
although  the  physics  of  its  formation  and  subse¬ 
quent  development  is  quite  different.  It  is  possible 
that  for  certain  gas  mixtures,  photoionization  may 
be  present;  however,  the  electron  kinetic  effects 
we  have  described  still  play  the  fundamental  role. 

Once  the  avalanche  front  makes  contact  with  the 
anode,  the  fteld  in  the  cathode  side  of  the  avalanche 
is  greatly  enhanced.  This  happens  because  a  con¬ 
ducting  path  now  exists  between  the  cathode  tip  of 
the  avalanche  and  the  anode.  Moreover,  the  field 
lines  converge  toward  this  tip,  so  that  subsequent 
avalanches  emerging  from  the  cathode  can  com¬ 
plete  the  final  bridging  of  the  gap. 

As  the  applied  voltage  increases,  two  things 
happen:  The  avalanche  “bursts”  closer  to  the 
cathode,  and  the  number  of  untrapped  runaway 
electrons  increases.  These  runaways  have  been 
experimentally  observed  for  voltages  two  to  three 
times  the  self-breakdown  voltage. 

in.  MATHEMATICAL  FOUNDATIONS  OF  THE  TWO-GROUF 
MODEL 

From  the  discussion  in  Sec.  n  it  is  evident  that 
the  breakdown  of  a  gas  at  higher  overvoltages 
(i.e,,  V>20%  self-breakdown)  is  governed  by  elec¬ 
tron  kinetics  in  a  highly  inhomogeneous,  three 
dimensional  field.  A  full  treatment  of  the  problem 
would  require  a  knowledge  of  the  space- time-de 
pendent  electron  energy  distribution  function  and  of 
the  space- charge  fields  in  the  region  bounded  by 
the  electrodes,  subject  to  the  boundary  conditions 
defined  by  the  secondary  processes  occurring  at 
the  electrodes  and  by  the  external  supply  circuit 
This  problem  is  indeed  formidable,  and  the  only 
possibility  of  it  being  solved  is  via  computer 
analysis.  In  this  paper  we  will  not  give  either  an 
analytical  or  a  computational  description  of  the 
physics  presented  in  Sec.  H;  however,  we  will 
give  a  mathematical  justification  of  the  hypothesis 
presented. 

Our  point  of  departure  is  the  state  of  the  ava¬ 
lanche  just  prior  to  stage  IV  (see  Sec,  II  and  Fig. 

1).  At  this  point,  we  have  a  weakly  ionized  medium 
(t3rpicai  values  for  the  density  of  electrons  are 

particles/cm^)  in  the  presence  of  an  ex¬ 
ternal  electric  field.  From  Gurevich's  work^  we 
say  that  the  electron  energy  distribution  function 


can  be  thought  of  as  being  formed  of  two  groups: 
the  main  distribution  electrons  and  the  fast  elec¬ 
trons.  This  discontinuity  is  again  of  a  mathe¬ 
matical  nature  to  make  the  problem  tractable. 

Although  the  conditions  existing  during  break¬ 
down  are  transient  in  character  and  bounded  in 
space,  we  will  assume  that  the  fast  electron  dis¬ 
tribution  just  prior  to  stage  IV  is,  to  zeroth  order, 
given  by  the  quasistatic  distribution  obtained  by 
Gurevich^  for  the  case  of  a  weakly  ionized  plasma 
in  an  infinite  domain;  that  is, 

ff(v,  t)  =‘f{vo,  t)  exp^-^  jf  vdv  .  (3) 

where  Vq  is  the  lower  boundary  of  the  fast  electron 
distribution,  Vq  ~  (5e^/m)^^^;  c,  is  the  ionization 
energy;  F{v}  is  the  effective  retarding  force  on 
an  electron  moving  with  a  velocity  v  in  a  neutral 
gas;  zxidfiv^yt)  is  the  main  distribution  evaluated 
at  v  =  Vo  (and  it  may  be  taken  as  Maxwellian).  The 
fast  electron  distribution  if  valid  for  y>  Vg. 

Instead  of  the  analytical  expression  used  by 
Gurevich^  for  the  retarding  force  F(t;)  (see  Ref. 

21  for  further  discussion  of  this  retarding  force), 
we  retain  the  same  functional  dependence  on  v  and 
introduce  two  parameters  A  and  B  which  are  thea 
determined  from  experimental  data.^^  The  reason 
for  this  is  that  there  is  a  disagreement  between 
the  values  obtained  experimentally  and  analytically 
for  the  position  and  height  of  the  maximum  of  the 
retarding  force.  Thus,  we  take  F(v)  to  be 

F(r)=A>ln(t;VBOA/, 

or  letting u-i mv^  (4) 

F(?4)=:A/>  ln(M/S)/M. 

The  function  F(m)  is  sketched  in  Fig.  4.  We  can 
now  find  the  injection  energy,  Le.,  the  minimum 
value  of  energy  that  a  fast  electron  must  have  to 
run  away  from  the  avalanche  and  be  injected  into 
the  decreasing  field  region  ahead  of  the  avalanche. 
From  Eq.  (3)  the  runaway  condition  is  defined  as 

F(u)  =  eE  (5) 

and  yields  the  minimum  energy  for 
(see  Fig.  4),  where  is  the  maximum  field  at  the 
front  of  the  avalanche  (Fig.  2).  Electrons  with 
energy  can  run  away  from  the  avalanche. 
Gurevich^®  found  that  the  flux  of  runaways  falls  off 
rapidly  with  increasing  ratio  F^/eE^  where 

is  the  maximum  value  of  the  retarding  force, 
so  that  the  greatest  number  of  runaway  electrons 
will  be  found  near 

In  a  reference  frame  moving  with  the  main  ava¬ 
lanche,  the  energy  change  of  one  of  these  Injected 
electrons  per  unit  distance  away  from  the  ava¬ 
lanche  front  is  given  by 
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FIG.  4.  Effective  retarding  force  as  a  function  of 
electron  energy. 


du/dx- eE{x)  ,  (6) 

where  u  is  the  fast  electron  energy  in  eV.  At  .r 
=  0,  «  =  injection  energy  of  a  given  fast  electron. 

From  Eq.  (6)  and  from  Fig.  5,  we  note  that  an 
injected  particle  with  energy  will  start  to 
gain  energy  because  du/dx  | >  0.  However,  if 
the  rate  of  energy  gained  from  the  field  is  less 
than  the  rate  of  energy  loss  due  to  collisions,  the 
particle  will  slow  down  and  be  trapped.  This  is 
due  to  the  fact  that  the  particle  is  moving  in  a  re¬ 
gion  of  decreasing  field  with  a  larger  gradient 
than  that  of  the  loss  function  F. 

In  (tt,  x)  space,  the  equation 

eE{pc)^F{u)  (7) 

defines  the  locus  of  points  for  which  du/dx A 
particle,  starting  with  energy  Mq,  whose  trajectory 
crosses  this  line  will  become  trapped.  Thus  there 
is  a  range  of  injected  particle  energies  u^<u<u^ 
which  will  slow  down  and  become  trapped,  u^  is 
the  injection  energy  of  a  particle  whose  trajectory 
is  tangent  to  the  locus  curve  [Eq.  (7)]  for 
The  electrons  in  this  range  of  energy  become 
trapped  at  various  distances  from  ttie  avalanche. 
Thus  a  narrow  channel,  as  described  in  Sec.  II, 
evolves  from  the  main  avalanche.  To  show  this 
effect,  it  is  necessary  to  solve  Eq.  (6)  with  F(m) 
given  by  Eq.  (4).  To  do  so,  we  must  know  the  de¬ 
pendence  of  the  electric  field  with  x. 

It  is  common  practice  to  assume  that  the  avar 
lanche  may  be  replaced  by  either  a  spherical  cloud 
of  electrons^  or  by  a  conducting  sphere,^®  We, 
however,  will  assume  that  just  prior  to  stage 
IV,  the  avalanche  may  be  represented  by  a  con¬ 
ductor  with  an  excess  of  electrons.  The  reasoning 
is  as  follows:  electrons  after  traveling  a  cer¬ 

tain  distance  produce  ions  and  electrons.  The 
rti+nz  electrons  and  ions  act  collectively  to 
establish  a  minimum  field  configuration  inside 
the  avalanche.  This  configuration  can  be  thought 


of  as  being  formed  by  the  superposition  of  a  dipole 
charge  distribution  (i.e.,  similar  to  the  uncharged 
conducting  sphere  approach)  and  an  electron  layer 
(similar  to  the  electron- cloud  approach).  The 
axial  field  at  a  distance  x  from  the  boundary  of  the 
avalanche  may  then  be  written  as^^ 

where  is  the  radius  of  the  avalanche  just  prior 
to  the  transition  into  stage  TV  and  E^  is  the  applied 
electric  field.  As  the  applied  field  increases  (e.g. 
50%  overvoltages),  the  transition  to  stage  IV  oc¬ 
curs  closer  to  the  cathode.  As  this  happens,  the 
contribution  to  Eq.  (8),  due  to  the  (n^  - 1)  ions  left 
behind  by  the  avalanche,  must  be  taken  into  con¬ 
sideration.  Their  effect  is  to  reduce  the  value  of 
the  second  term  in  the  right-hand  side  of  Eq.  (8).  * 
At  higher  fields  (^80%  overvoltages),  the  transi¬ 
tion  occurs  very  close  to  the  cathode.  For  these 
cases,  the  whole  space  charge  has  a  dipole  struc¬ 
ture  so  that  the  second  term  in  Eq.  (8)  disappears. 

Numerical  solutions  of  Eqs,  (4),  (6),  and  (8)  are 
shown  in  Fig.  5  for  the  case  of  breakdown  in  a  1- 
cm  gap  in  nitrogen  at  atmospheric  pressure.  The 
dotted  curve  is  the  locus  curve,  i.e.,  Eq.  (7).  The 
differences  between  Figs.  5(a),  5(b),  5(c),  and 
5(d)  are  in  the  values  of  and  Eq  used.  In  Fig.  6, 
we  show  tiie  minimum  injection  energy  as  a 
function  of  the  maximum  electric  field  for  this 
example.  The  values  for  the  parameters  A  and  B 
in  Eq.  (4)  were  found  from  the  data  in.®^ 

The  results  are  qualitatively  similar:  Injected 
fast  electrons  with  energies  in  the  range  u^<u<ut 
become  trapped  at  various  distances  ahead  of  the 
avalanche.  Electrons  with  u>Uf  will  actually  run 
away.  Their  number  increases  as  the  applied  field 
increases.  When  the  runaway  flux  is 

substantial.  Under  these  conditions,  the  fast  elec¬ 
tron  distribution  is  no  longer  given  by  Eq.  (3).  A 
solution  to  the  Boltzmann  equation  when  strong 
electric  fields  are  present  must  be  obtained. 

The  minimum  injection  energy  and  the  trapping 
distance  are  a  function  of  and  Eq,  The  conse¬ 
quence  of  this  dependence  is  discussed  in  the  next 
section. 

IV.  DISCUSSION 

It  is  impossible  to  apply  a  fundamental  model 
alone,  as  presented  in  Sec.  Ill,  to  a  specific  ex¬ 
periment.  In  general,  each  experiment  involves 
the  fundamental  physical  processes  and  processes 
related  to  idiosyncrasies  of  a  particular  experi¬ 
ment  It  is  sometimes  difficult  to  isolate  the  con¬ 
tribution  from  each  process  thus  leading  to  mis¬ 
interpretation.  In  the  early  days  of  the  streamer 
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FIG.  5.  EvolutiOQ  of  injected  electron  energy 
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FIG.  6.  Minimum  Injection  energy  as  a  function  of 
maximum  electric  field. 


model,  Dickey^^  questioned  the  necessity  for  in¬ 
troducing  a  “new”  model,  besides  the  Townsend 
model,  to  explain  the  short  time  lag  experimentally 
observed  by  Rogowski“  and  later  by  Fletcher.^*^ 

He  argued  that  these  time  lags  were  due  to  the 
particular  experimental  arrangement  and  not  due 
to  any  new  physics  in  the  breakdown  process. 

In  the  discussion  of  our  model,  we  will  con¬ 
centrate  on  the  general  structure  of  the  observed 
breakdowns.  We  will  also  examine  the  transitions 
that  may  occur  in  any  specific  experimental  ar¬ 
rangement  due  to  the  fundamental  processes  in  the 
gas  alone. 

In  the  regime  of  interest  to  our  model  (i.e., 
^>20%  self-breakdown),  two  kinds  of  breakdown 
development  have  been  observed:  (1)  filamentary 
breakdown  and  (2)  broad  breakdown.  Broad  break¬ 
down  has  been  observed  either  at  very  large 
voltages^"^  iV>2  times  self -breakdown)  or  at  lower 
voltages  when  the  initiating  electrons  are  generated 
over  a  large  surface  of  the  cathode.^^  In  all  other 
cases,  filamentary  breakdown  occurs. 

The  principal  parameters  that  we  need  to  de¬ 
scribe  the  structure  of  the  breakdown  are  (1)  the 
applied  electric  field  (2)  the  spatial  dependence 
of  the  electric  field  just  ahead  of  the  avalanche, 
and  (3)  the  frictional  force,  F^k),  on  an  electron  as 
a  function  of  its  energy  for  a  given  gas.  Note  that 
F(u)  also  depends  on  the  pressure.  From  Fig.  5 
we  note  that  if  the  avalanche  is  initially  small  in 
cross  section  (this  occurs  when  the  electrons 
initiating  the  avalanche  are  generated  in  a  small 
area  of  the  cathode),  the  inhomogeneity  in  the 
electric  field  is  large.  This  implies  that  the  ener¬ 
gy  spectrum  of  the  trapped  electrons  is  broad, 
i.e,,  many  electrons  will  be  trapped.  Moreover, 
these  electrons  will  be  injected  along  the  axis  due 
to  the  “escape  cone”  phenomenon.  This  behavior 
corresponds  to  the  range  of  parameters  where 
filamentary  breakdown  has  been  observed. 


As  the  field  is  increased,  we  reach  the  regime 
described  by  Mesyats'  avalanche  chain  model.^® 

He  proposed  that  the  main  avalanche  slows  down 
and  that  a  new  avalanche  develops  from  the  pre¬ 
vious  one.  This  process  is  repeated  a  number  of 
times  until  an  avalanche  chain  is  formed. 

In  our  model,  the  process  is  continuous.  At  the 
higher  field  values,  the  width  of  the  energy  spec¬ 
trum  of  the  trapped  electrons  decreases  due  to  the 
fact  that  the  decrease  in  the  injection  energy 
will  be  smaller  than  the  decrease  on  the  lower  en¬ 
ergy  limit  for  runaway,  i.e.,  u^.  Accordingly,  the 
range  of  the  trapped  electrons  decreases.  The 
avalanche  is  extended  via  these  short-range  elec¬ 
trons  [see  Fig.  5(a)].  The  number  of  runaway 
electrons  and  their  effect  in  this  regime  is  still 
small.  However,  as  the  field  increases  further 
(Le;^2  times  self-breakdown),  the  injected  elec¬ 
trons  are  no  longer  trapped.  In  this  regime  eFCr) 
>F(u{x))  and  electrons  continuously  gain  energy; 
they  all  run  away  [compare  Figs.  5(a)  and  5(d)]. 
This  is  the  regime  discussed  by  Babic  ef  Our 
model  is,  in  this  regime,  similar  to  theirs;  how¬ 
ever,  we  do  not  assume  that  the  injected  electrons 
have  zero  initial  velocity  and  we  are  consistent  in 
our  definition  of  runaway  electrons.  This  dis¬ 
charge  in  this  regime  has  a  diffuse  or  multi¬ 
channel  character. 

A  situation  similar  to  the  above,  i.e.,  no  trapped 
electrons,  exists  when  the  breakdown  is  initiated 
by  a  large  number  of  electrons  distributed  over  a 
large  area.^'*  In  this  case,  the  inhomogeneity 
ahead  of  the  avalanche  is  small  so  that  once  again 
injected  electrons  actually  runaway  [see  Fig.  5(c)]. 
The  result,  as  in  the  very  high  field  conditions,  is 
a  broad  breakdown  due  to  ionization  and  radiation 
generated  by  these  runaways, 

A  subject  of  interest  has  been  whether  or  not  an 
abrupt  decrease  in  the  formative  time  would  be 
observed  as  the  percent  overvoltage  above  self¬ 
breakdown  is  increased.  In  this  respect  the  re¬ 
sults  have  been  mixed.^®  In  some  cases,  in  oxy¬ 
gen,  and  nitrogen,  for  example,  a  decrease  was 
observed  for  large  gap  spacing.  In  hydrogen,  how¬ 
ever,  it  was  not.  This  behavior  can  be  explained 
by  considering  the  magnitude  of  the  maximum  of 
the  retarding  force  for  each  gas  and  the  percent 
overvoltage  used.  The  number  of  injected  elec¬ 
trons  is  a  function  of  where  (see 

Fig.  4),^®  If  (see  Sec.  Ill)  the  number  of 

injected  electrons  is  practically  zero  so  that  no 
transition  will  be  observed.  As  E^  increases,  the 
number  of  injected  electrons  increases.  Accord¬ 
ingly  a  continuous  change  in  the  slope  of  the  time 
lag  versus  percent  overvoltage  curve  may  be  ob¬ 
served.  This  was  observed  experimentally  by 
Phillip  and  Allen.^®  We  were  not  able  to  check  our 
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model  for  all  the  gases  they  used  due  to  lack  of 
data  on  the  frictional  force.  For  the  cases  of  hy¬ 
drogen,  oxygen,  and  nitrogen,  for  which  data  is 
available,  the  observed  time  lag  versus  over¬ 
voltage  curve  may  be  explained  by  our  model 

In  the  case  of  hydrogen,  the  applied  overvoltages 
used  were  too  low,  so  that  no  change  in  the  slope 
was  observed*  Tlie  maximum  of  the  retarding 
force  for  hydrogen  is  relatively  high,  while  the 
breakdown  voltage  for  a  1 -cm  gap  is  small  (i.e., 
*15  kV).  Therefore,  large  percentage  overvol¬ 
tages  (“3  times  self-breakdown)  will  be  required 
before  injection  can  begin. 

For  the  cases  of  oxygen  and  nitrogen,  changes 
were  observed  for  a  3- cm  gap  at  atmospheric 
pressure.  The  changes  were  more  pronounced  in 
oxygen  than  in  nitrogen.  The  maximum  retarding 
force  for  oxygen  is  less  than  for  nitrogen;  more¬ 
over,  the  breakdown  voltage  for  oxygen  is  higher. 
This  combination  may  account  for  the  observed 
differences. 

As  the  gap  width  decreased  for  the  same  range  of 
overvoltages,  the  change  in  the  slope  decreased. 
For  a  gap  width  of  1  cm,  no  change  was  observed 
In  nitrogen  for  up  to  25%  overvoltages.  As  for 
hydrogen,  these  values  are  again  too  low  for  in¬ 
jection  to  occur.  No  experiments  were  performed 
for  oxygen  in  1-cm  gaps  although  the  trend  seems 
to  indicate  that  a  similar  behavior  would  be  ob¬ 
tained.  The  changes  observed  for  3-cm  gaps  in 
these  two  gases  at  these  lower  overvoltages  may 
be  due  to-  large  space-charge  distortion  which 
lower-s  the  injection  energy.  This  may  occur  if 
the  magnitude  of  the  ‘tree  charge”  just  prior  to 


stage  IV  is  large  (see  Sec.  Ill),  in  which  case  the 
second  term  in  Eq.  (8)  will  be  larger. 

The  next  step  in  the  development  of  the  model  is 
to  do  a  numerical  analysis  of  a  fluid  model  for 
various  particle  species.  The  important  point  is 
that  the  electron  fluid  must  be  divided  into  two 
groups:  the  fast  and  slow  electrons,  A  number  of 
experiments  should  be  carried  out  to  test  the 
validity  of  the  theory  and  its  basic  assumptions. 

In  particular,  since  the  retarding  force  in  Eq,  (6) 
is  a  function  of  the  properties  of  the  gas,^^  the 
appearance  of  the  runaway  electrons  will  depend 
strongly  on  the  gas  used.  A  series  of  experiments 
should  be  carried  out  to  test  this  dependence. 
Moreover,  experimental  verificatioa  of  the  time 
of  appearance  of  the  runaway  electrons  as  a  func¬ 
tion  of  percent  overvoltage  should  be  provided. 

ACKNOWIEDGMENTS 

We  wish  to  express  appreciation  to  Professor 
G.  A.  Mesyats  and  Professor  M.  O.  Hagler  for 
stimulating  discussions,  and  to  Dr.  B.  D.  Black- 
well  for  his  assistance  with  the  computer  graphics. 
We  also  want  to  thank  the  reviewer  for  suggesting, 
among  other  things,  changes  in  the  title  of  the 
paper  to  its  present  form.  This  work  was  support¬ 
ed  by  the  Naval  Surface  Weapon  Center  under  Con¬ 
tract  No.  N60921-79C-A187  and  in  part  by  the  Air 
Force  Office  Scientific  Research  under  Contract 
No,  AFOSR  76-3124A.  The  findings  in  this  paper 
are  not  to  be  construed  as  an  official  Department 
of  Navy  position  unless  so  designated  by  other 
authorized  documents. 


Llewellyn- Jones,  Ionization  and  Breakdown  in  Gases 
(Methuen,  Loodon,  1966),  p.  30. 

Raether,  Electron  Avalanches  and  Breakdown  in 
Gases  (Butterworths,  London,  1964). 

^J.  Dutton  and  W.  T.  Morris,  Br.  J.  Appl.  Phys.  !£,  1115 
(1967). 

^L,  H;  Fisher  and  B.  Bederson,  Phys.  Rev.  109 
(1951). 

®A.  J.  Davies,  C.  J.  Evans,  and  P.  M.  Woodison, 

Comput.  Phys.  Common.  287  (1978). 

®A.  J.  Davies,  C.  J.  Evans,  and  F.  Llewellyn- Jones, 
Proc.  R,  Soc.  London  A281,  164  (1964). 

^A.  L.  Ward,  Phys.  Rev.  1357  (1965). 

*G.  M.  Petropouios,  Phys.  Rev.  78,  250  (1950). 

®E.  D.  Lozanskii,  Zh.  Tekh.  Fiz.  1563  (1968)  [Sov. 
Phys.— Tech.  Phys,  1269  (1969)1. 

Rogowski,  Archiv  Electrotech.  (Berlin)  20’ 

(1928). 

C.  Fletcher,  Phys.  Rev.  7^,  1501  (1949). 

^^F,  Llewellyn- Jones,  Ionization  and  Breakdown  in  Gases 
(Methuen,  London,  1966),  p.  149. 


a.  Mesyats,  Y.  I.  Bychkov,  and  .A.  1.  Iskordskii, 
Zh,  Tekh.  Fiz.  38,  1281  (1968)  [Sov.  Phys.— Tech. 
Phys.  13,  1051  (1969)]. 

^^F.  R.  Dickey,  Jr.,  J.  Appl.  Phys.  23,  1336  (1932). 

^*Y.  L.  Stankevich  and  V.  G.  Kalinin,  Dokl.  Akad.  Nauk 
SSSR  177,  72  (1967)  [gov.  Phys.— Dokl.  1042 
(1968)!. 

^*L.  P.  Babich  and  Y.  L.  Stankevich,  Zh.  Tek.  Fiz. 

1669  (1972)  [Sov.  Phys.— Tech.  Phys.  1£,  1333  (1973)]. 
^^L,  V.  Tarasova,  L.  N.  Khudyakova,  T.  V.  Loiko,  and 
V.  A.  Tsukerman,  Zh.  Tekh.  Fiz.  U,  564  (1974)  [Sov. 
Phys.— Tech.  Phys.  351  (1974)]. 

^*V.  V.  Kremnev  and  G.  A.  Mesyats,  Zh.  Prikl.  Mekh. 
Tekh.  Fiz.  1  [J.  Appl.  Mech.  Tech.  Phys.  (USSR)  40 
(1971)1. 

^*E.  D.  Lozanskii,  Usp.  Fiz.  Nauk.  117,  493  (1975)  [Sov. 
Phys.— Usp.  1^,  393  (1976)]. 

V.  Gurevich,  Zh.  Eksp.  Teor.  Fiz.  1296  (1960) 
[Sov.  Phys.— JETP  12,  904  (1961)1. 

D.  Landau  and  E.  M.  LifShitz,  Electrodynamics  of 
Contimious  ytedia  (Pergamon,  London,  1966),  p.  344. 


42 


DEVELOPMENT  OF  OVERVOLTAGE  BREAKDOWN  AT... 

”L.  R.  Peterson  and  A.  E.  S.  Green,  J.  Phys.  B  1131  Koppitz,  J.  Phys.  D  6,  1494  (1973). 

(1968).  “  25k,  R.  Allen  and  K.  PUlUps,  Proc.  Roy.  Soo.  London 

^J.  A.  Stratton,  Electromagnetic  Theory  (McGraw-Hill,  A278.  188  (1963). 

New  York,  1941),  p.  205. 


43 


THIRTY-THIRD  GASEOUS  ELECTRONICS  CONFERENCE 


PLEASE  TYPE  NAME  &  ADDRESS 

E.E.  Kunhardt _ 

Iflillzed  Gas  LaboratorY _ 

Texas  Tech  University 

P.n.  Rnx  4439 _ 

Lubbock.  Texas _ 79403 _ 


7-10  October  1980 
Norman,  Oklahoma 

DO  NOT  WRITE  IN  THIS  SPACE 

Serial  No - 

Accepted:  Yes  - No  - 

Session _ 

Number - 

Date  Conf.  - 


AUTHORS  PLEASE  NOTE 

Indicate  in  what  subject  category  your 
paper  best  belongs  (see  preliminary 
announcement  for  category) 

1  2  3  4  5  6 

8  9  10  11  12 

13  14  15 


Development  of  Overvoltage  Breakdown  at  High 
Gas  Pressure  -  E.E.  KUNHARDT  and  W.W.  BYSZEWSKI,  Texas 
Tech  University*  --  A  model  for  the  development  of  elec- 
trical  breakdown  in  high  pressure  gases  is  presented.  It 
describes  the  initial  phases  of  breakdown  in  the  regime 
where  the  Townsend  avalanche  mechanism  does  not  apply. 

In  the  model,  the  energy  distribution  function  for  elec¬ 
trons  in  the  advancing  avalanche  is  assumed  to  have  two 
components:  The  fast  and  the  slov/  electrons.  The  fast 
electrons  “runaway"  from  the  avalanche,  but  are  "trapped" 
at  various  distances  ahead  of  the  avalanche,  depending  on 
the  energy  of  the  electrons  and  the  properties  of  the 
avalanche  and  the  background  gas.  This  mechanism  plays 
the  fundamental  role  in  increasing  the  speed  of  propaga¬ 
tion  of  the  avalanche  when  the  applied  voltage  is  ^  20 
percent  of  the  self  breakdown  voltage.  In  light  of  this 
model  a  brief  discussion  of  experimental  results  if 
given. 
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Numerical  Solution  of  Poisson's  Equation 
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Objective:  to  solve  :-p  3^  (r  3^  4>)  +  3^  =  -4irp  (1) 

with  boundary  conditions: 


(j)  =  0 

0 

z  =  0 

(cathode) 

(j)  =  0 

@ 

z  =  d 

(anode) 

3^  (j,  =  0 

0 

r  =  0 

(^  =  0 

0 

Where  r^^  is  the  "furthest  radius"  of  interest.  (2) 

Motivation:.  The  above  "mathematical"  problem  is  applicable  to  the 
"physical"  problem  of  finding  the  values  of  electric., 
field  in  a  gap  due  to  the  space  charge  in  the  gap.  Cyl¬ 
indrical  symmetry  (which  is  not  a  major  restriction)  is 
assumed. 

Importance:  A  major  obstacle  in  tracking  the  space-time  development 

of  the  breakdown  of  an  overvolted  gap  is  the  time  (com¬ 
putational)  required  to  solve  the  Poisson  equation,  A 
fast  alogarithm  will  result  in  large  savings  of  time, 
which  in  turn  implies  that  we  will  be  able  to  track  the 
development  of  the  breakdown  for  a  longer  period  of  time 


Results:  We  have  developed  a  direct  alogarithm  for  solving  Poisson's 


Equation  in  Cylindrical  coordinates. 

Technique: 

1.  Apply  Fourier  transform  in  z,  making  use  of  the  FFT  alogarithm. 

2.  Use  Cubic  Splines  to  solve  for  the  fourier  coefficients  as  a 
function  of  r. 

3.  Directly  obtain  the  electric  field  by  inversion  of  appropriate 
spline  coefficient. 

Description  of  Technique: 

Let 

<|)(r,z)  =  ^  (i>  ,^(r)  sin  ^  z 
and 

p(r,z)  =  ^  P|^(r)  sin  ^  z 

Substituting  into  Equation  (.1)  and  using  the  orthogonality  properties 
of  the  sine  functions  in  the  interval  0  £  z  _<  d,  we  get: 

4  \(r)  +  1  ^  ♦,(r)  -  (  )2  4,^Cr)  =  -4.p,(r)  (3) 

dr 

Equation  (3)  is  a  set  of  k  equations  since  we  only  need  k  harmonics  i.e. 
equal  to  the  number  of  grid  points  in  z. 

The  boundary  conditions  that  go  with  Equation  C3)  are: 

3^  (j),^  =  0  @  r  =  0 

Since  the  developing  avalanche  has  a  very  small  radius  .1  cm),  it  is 
desirable  to  expand  the  region  around  the  origin  (r  =  0)  by  using  a  trans 


formation  of  coordinates  in  the  r  direction.  A  transformation  which  will 
suit  our  purpose  is: 

?  =  a{l  -  e’^*")  (4) 

Where  a  and  b  are  parameters  chosen  so  as  to  provide  the  desirable  ex¬ 
pansion  around  the  axis.  Note  that  for  o  £  r  £  <»;  the  r,  range  is  0  £  c  £  a- 
Using  the  transformation  in  (4),  Equation  (3)  becomes: 


^  (f).  U)  +  p(c) 
dC 


_d 

dc 


\(?)  -  q(?)  =  -s(?) 


Where 


P 


In  (  ^)  +  1 

(?)  = - 

In  (5fi) 


1 

a-Z 


(5) 


q(?)  =  +  -T - ^2  and  SU)  = 

bHa-zr 

Equation  (5)  is  valid  everywhere  except  for  z  = 
?  ->■  0,  we  obtain: 


4TrPk 

(a-z)^ 

0.  Taking  the  limit  as 


2b  V, 


-27rPk(0) 


h2,2 

b  a 


(6) 


To  solve  Eqs.  (5)  and  (6),  define  the  cubic  spline  as  a  function  G  (?)  s  C2 

r2l 

of  the  form,  in  each  of  the  subintervals  [?^. ,  ■* 

where  G„  .  are  constants. 
n,i 

Now  we  define  the  system  of  Eqs.  for  G  . ,  n  =  0,  1,2,3;  i  =  0,1, - N; 

n  j  1 

such  that  Eq.  (7)  satisfies  Eq.  (5)  and  (6)  at  the  points  . . 

and  G(?),  G'  (?),  G"  (?)  are  continuous  at  the  points  z^, . . 

From  Eq.  (7),  we  get: 


(«1)  =  '=0.1 


4 


(8) 

(8b) 

(8c) 


The  grid  structure  in  ?  is  shown  below. 


Si_i  s. 


?1  ?2 


i-1  1 


Applying  continuity  requirements  at  ith  point: 

il  S.(s.j)  =  (9) 

Noting  that:  o 

i-1 

^i-l^)  =  ^o.i-1  "  '^l.i-l  ^  ^S-.l  ^  ^2.1-1  -T-^  S,i-1  (lOa) 


^?'^i-l^*’i^  ^  ‘^l,i-l  *^2,i-l  ^2, i-1  ^^i-1  S,i-1  ■ 

where  A?^_i  =  -  ?.j_i 

Using  (8a),  (9),  and  (10a),  ^  ^^3 

%.i  =  ^0,1-1  ^  *  ^2, i-1  ^  "  ®3,i-l  ^ 

Using  (8b),  (10b),  and  (12): 

2 

^1,1  ='^1,1-1  *®21-1  ‘'1-1  *®3.1-1 


(10b) 


c)_  9^^  6.(;.)  =  9^^  (?^.)  (14) 

Using  (8c),  (10c),  and  (14) 

®2,1  '  '^2,1-1  ^  S,i-1  ‘'i-l 

Equations  (11),  (13),  and  (15)  are  valid  for  the  points  i  =  1  •— N. 
Besides  these  continuity  equations,  the  spline  functions  must  satisfy 
the  differential  equation  at  each  point.  From  equations  (5)  and  (6)  and 
equations  (8): 
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(16) 


^2,i  ‘^1,1  ^  ^  ° 

The  above  equation  is  re-written  as: 

'^2,i  ^  Pi  ®l,i  ■  Pi  ®o,i  "'^• 
where  =  p{c^)  i  f  0 

=  0  i  =  0 

q^-  =  q(?i-)  i  =  0,  ...  N  +  1 

s.  =  s(c^.)  i  =  0,  ...  N  +  1 

Using  Equations  (11),  (13),  (15),  and  (16),  we  obtain  the  following  Equation 

for  G-|  .  (after  considerable  algebra): 

*  »  • 


A,  ]  G, 

,i+l 

II 

1 — 1 

m 

i+l  ^i^i 

i+l 

Bib.  + 

*-i+l 

C,d, 

^  '^l,i-l 

i+1 

C,A,  - 

^i+1 

i+l 

E .  b .  + 
11 

^i+1 

C,fi 

1 _ 1 

CO 

1 

for 

i  =  1, 

•  •  #  • 

•  >  N 

(17) 


Note  that  from  boundary  conditions  and  equation  (8b): 

,0  ^  ^1  ,N+1  ^  ° 

That  is,  the  radial  field  is  zero  on  the  axis  and  at  5^^. 

Equation  (17)  is  a  tridi agonal  system  of  equations  for  G-j  ^  which  can  be 
inverted  using  standard  computational  techniques.  The  constants  A,  B,  C, 


D,  E,  and  a  ,  b,  c,  d,  e,  f  are  given  by: 

A?i.i  q^C. 


A.  =  1  + 


^^i-1 


Pi  - 


B.  =  1  - 


A? 


i-1 


AS'-l 


^^.d. 

1 


'i-1 


A?.; 


'^i  = 


= 


1 


2 

i-1 


Pi-1  ^ 


^?i-l 


Pih 


i-1  Pi®i 


^i  = 


A?i.i  ACi.i  Pi^- 


a.  =  1  + 
1  1 


b.  =  1  - 


i-1 

6 

2 

i-1 


^i 


fi-i 


*^1  ^ 


i-1 
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d..  = 


A?i.l 


'i-1 


f..  = 


i-1 


6 

-1 

^  1-1 


Note  from  equation  (8b)  that  G,  is  the  fourier  transform  of  the 

radial  Electric  field,  so  we  need  not  find  any  of  the  other  coefficients. 
Inverting  G,  ^  (?,•),  we  obtain  the  radial  field.  The  inverse  cosine  trans- 

I  j  1  I 

form  of  K  G  •  gives  us  the  axial  field. 

0  j  I 
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Sine  and  Cosine  transform  subroutines 
and  Spline  subroutines. 


The  subroutines  FFTCN,  X,  Itypel,  Itype2)  and  spline  are  used  to  solve 
the  Poisson  equation  in  cylindrical  cooordi nates  with  axial  symmetry.  The 
procedure  is  as  follows: 

Given  a  space  charge  p(r);  the  forward  sine  transform  is  taken  (i.e, 
Itypel  =  1,  ItypeZ  =  -1).  Then  the  radial  equation  (Eq.  3  of  outline) 
is  solved  for  each  fourier  component  using  subroutine  spline.  The  space 
charge  is  brought  into  the  subroutine  as  ER(I,  K).  Use  is  made  again  of 
the  FFT  to  obtain  the  radial  field  and  the  longitudinal  field  at  each  radial 
point,  as  a  function  of  z  (axial  coordinates),  via  inverse  transformation 
(Itype2  =1).  Note:  The  Spline  subroutine  shown  does  not  use  coordinate 
transformation.  For  this  case  p^C?)  in  Eq.  5  is  given  by  1/r^.  and  q^(c)  is 
given  by  (-^)^ 
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i 

i 

i 

)i 

'V _ c 
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SUBSOUTIME  FFT (N»X* TTYHEl « TTYPF?) 

r 

subroutime  for  ftt  used  iNTn:wi^rrTN(r3TNEn3R  corsT!vrE“TiR7m!SF0TrMs 

M:;  C 

note  :  IE.  signal  HAS  (2**(N-1)*1.  POINTS 

i.:  C 

H . c- 

. C 

TTYPEi! . -1  FOP  COS  TPANSFORN*  i  FOP  SIN  TRANSFORM . . . . 

U:  ■  -  C 

Isr'V',  C 

U  '  ■  c 

i  ■  •  - 

TTYPEH'I  F'FOR"' inverse  TRANSFORMf  -1  FOR  FORWARD — MRA NSFORK - ^ - 

TVRE  XOT^PrETK  W^»FBUF*  ZTC^^UXT7^ 

COMMON/FT/W (?56) »W1 (256)  . 

OTMENSTOM  FRUF (512) ,2(512) ♦X (IPR)  . .  . 

!:/-  ANsNT 

h  NH=NT/?  :  :  ' 

c 

NKIrNH+l 

C 

!  '■ 

CHECK  IF  FORWARD  OR  INVERSE 

!,  ■ 

10 

IFdTYPEp.GE.O)  10,40 

continue 

jZ... 

c 

20 

IFTT TY  P  ET“  .  G  F.  0 )  20',  A-O" 

CONTINUE . . . . .  . . 

c 

c 

TPANSFER~D“AT7r  TCT  F80F~Cr0FiPLEX  ARRaV 

FBUF{1)SCMPLX(0.,X(1) ) 

DO  30 

F3UF(T)=CWPLX(0.»Xa)  ) .  . . 

FRUF(NT+B-T)=-F8'JF(I)  . . 

30 

40 

CONTI  we 

GO  TO  80 

CONTINUE 

FRurci ) sc^'PLXrX (I )  ,(5.r  -  .  -  - 

DO  no  I=2»MH1  . 

FRlJF(n=C‘^°LX(X(I)  ,0.)  . . 

r 

c 

c 

CHECK  IF  SINE  OR  COSINE 

390 

IF ( I TYPEl • GE , 0 )  390,400 

COMTTMUE 

FPUF(mT*?-I)s-FBUF(I)  . 

c 

c 

c 

NOTE  THAT  FBUF (NT/I+1 ) =0 . 

GO  TO  110 

400  FBIIF  (MT+p-I  )  sFBUF  a  ) 
IIU  COMTIMIJF 

B  O  COTMTI  MOF - -  - 

C 

C  UO  STUFF 


on  ;35n  k=i,n 


FTNfS,5 


f  K-1  ) 

L  =  ?-»M 


11/12/80 
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C 

~  Trrr-?-r5~i^T7WT,T:“- - - 

iLM=r/L*M 
C  ' 

OTTTO'O  ^ 

. 2  0  0  Z  (  I  +  J- 1  j  =  F  fid  r  { I L  W  t.  J )  ♦  FB  U  F  ( TC  M  +  J  H )  . 

C . . . . 

■  2:it;:a=i»w 

■  c  -  ■■■.. 

~C  check  ’  t'P  '  FOR'^A'PO  OR  inverse -  - 

"C . . . . . . . . . . . 

. IF(  ITyPE? . .'GE. . 0) . 4|o-,-420  ’~ . 

4'2'0 Zn+J+M-i  Is  (F8IJF  ( lUHf  J 1  -FBUF  ( Il^M* )  »W  ( ILW+l)  ' 

60  TO  210 

410  Z(  l  +  J  +  M-l  )-(FedF(  ILH+J)  -FRliFf  lLM  +  J  +  NM)  J<»W1  (ILM+I  ) 
Z'l’lT’COKlT’TTruE  ^  ~  ' 

215  CONitTNiJE . . . . . . 

DO  220 . T=lVMt . . . . 

T2Q  FSUFfT  1=2(1) 

250  CONTTNUE 
C 

“C  ^CRETCTrTr'-FCf^prffR  D“0P"TNVEnT5E — ^ ^ - 

■  0 . ; . . . . . . . . . . 

00  600  t  =  l'*NHl . . . . . . . . 

- IFTTTYPEr27GETOT?T(n:2FO - - - ^ - - 

270  CONTINUE 

F8UF(I)sFBUF(r)/AM 
GO'~ro''54n  . 

280  CONTINUE . . . . . . . . 

r  . . . . . . . 

’X - CNECK“'rF~STNE ’"0 R  XO STNF - ^ - — - - 

C 

IF( ITYPEI .6E.0) 530f540 

. ~‘"530  X  rTTsxrMAOTFcrijFiryi  ~ 

GO ' TO  600  . 

540  X(T)=RFAL(FflllF(I  J  ) .  . 

6rrO~XOMTT.NnE  - - 

RETURN 

END 


c 

c 


S'JR'POUTTNiE  'SPLINE 


FTNi5.5 


11/12/90 
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C 

C 

C_ 

"C 

c 


c 

c_ 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 


c 


-  SUP90UT TME  EV ALUATES  THE  PAQIAL  EQUATION  USING  SP LINES 

ALL  CONSTANTS  NEEDED  IN  THIS  SURROUTINE  ARE  CALCULATED 
JM  “fHE  MAIN  PR06fiAM“Af'TirAR'E  COMMW“WrnH“Tl9lT''SWH 


D I M E N S  ION  P H I  ( 64  »  64 )  »  E R  ( 64  »  6 4 )  »  E  Z  { 6 4  ♦  6 4 ) . . . 

DIMENSION  ORGl  ( 128)  »0Re2  ( 1 28  )  » 0R^r?8T“ 

DIMENSION  alpha (128) *SI ( 128) »B M28) 

EQUrVALENCECPHIfEZ)  • 

COMMON  ERtE'Z  ^ - 

•■CnMMON/A/PO»ZO«NR»NZ/e7DRVDR7»ORSiDRS6 . . . . . . . - . . . 

. C  O  M  M  0  N  7  c  /  N  P  t ,  N  P  M 1 1 N  R  M  2 .  . . . . . . . . 

COMmon/o/  ITV*  PT"  ■ "  ;  “  ■  . — ^ 7: - — - 

P0  IS  THE  LENOTH  IN  THE  «  DIRECT ION»  ZO  IS  THE  LENGTH  INTHE 
"inoiRECTlON,  ^ 

'DP . t  g-'-fHE . STEP . SIZE  ■■"IN . R . - . - . - . 

DR2sD'R72",' . . . . . . . . . . 

DRSsDR^DR  - - - - - - : - : - — — 

'.ORSA*ORs/6r 

NR+I  -IS:  the  total  NUMBER  OF  POINTS  IN  R, 

MRlsMR  +  1  - — 

. NpMi=NR-'l: . . . . . . . . . . . . . 

N  R  M  PsN  R  -  2  - . .  . . - . .  . . . . . - . . . .  - . 

ITV  ’  1 5  THE  ETf5EN^  'VALUE  IN  ThTTZ"  0 1 RECrrDfrTnET'T^TN~smEnp'PA^ 


EVaITV*PT/ZO 


D  O  L  a  S  T  p  0 i  n  T I N  T  R I OT  A  G  0 N A  L  FACTOR  I Z A  T I O  N 


laAlRMl 

ORST0=(I-l )»D» 

'■ORSTi-^AOR - 

DPSTs ( r+t ) ^DO 

P?=l ./orst  . 

TTsTTZOPSTT - — ■ 

po=i,/n.PsTo 

Q1=EV«-EV 

“OTTsOl— - 

02=0  T 

■AST=  r;  -OTtfDRSS - 

8S1=1  ,>0fi*2.»DRS6 

CS1=0RS6*P1 

D  SI  =0  D‘^2  .»■  0  R"?  6  *  PO 

Esr=nps6 . 

FSlsp.^OPSS 


Al.lsl  .+D02*®!  ■fCSl*01*OR2/ASl 
BL]*1  .-0PP*Pn4.0Sl»DR2*Ql/ASl 
C  L 1  =  OP  P’*  0  0  9  S 1 »  0  R  2  «  0 1  /  A  S 1 
nLl=nP2+ESl*DR«*QI/ASl 


FTNsS.5  ■  ^  ■  . .  . .  ‘  11/12/80 


EL1=0P?+FS1*0R2*Q1/AS1 


s _ 

C 

AS2=1 .-Q?»nPS6 

k/ 

SS2  =  1 .♦2,*nPS8*01  "  ■ 

CS2=nRS6»P2 

DS2=0R-2.»0RS6*p1 

C 

'ES'2=ES1  .  '  ■  ■ 

ES?=FSI  . . . . 

k 

ffL?=l.*DP?#P2>CS2'»W2^‘277f5'2  ^ . . . . . . . 

eL2=l .-0R2*Pl*nS2»DR2*Q2/AS2 

Ct.2snR2»0l.  ♦BS2«DR2«Q2/AS2 

c 

0l2sOR?*ps2*Dp^<*(327AS2 

EL2=0R2-t-FS2*nR?*Q2/AS2 . . . . . .  . 

i'j:. 

■  A=ATl'Wir27Ct'2^Sr«'Ain/'CCT='^ - - - - 

8 (N) sOSl-3Sl*RLl/CLI 

5^0"  sF  <n  -  8  S 1 »  6X 1  /  C  L 1 

SI  =  ASl»EL2/CL2*ESl-BSl«r)Ll/CLl  . 

F2=Sf»ER  (NR,  TTV)  ♦S0*EH  (MRMl  f  ITV) . - . 

1- . ' 

~^PBA  (MysA  . 

.  SI'CNlsF?,.  .. 

;  GAMMA:50„:'  -  ■  : 

2r! 

ORBl  (N+l )  sAL2/C'L2 

. 0 R'B  2  ( N  + 1")  =  -B u 27 C L  2 . . . . . . . . 

0RB3  ( N  +  1)  =:EL2«'ER  (NRMl,  I  TV )  /CL2 . . . 

■'  • 

c 

c 

c 

evaluatf  interm^  the  tridiagonal  equation 

DO  2o  1 0  ■t  =  2»NPH2 

N  F  s  R  -  J 

PI=pn . . . . . . - . 

QIsQO 

asrsasi .  ■ 

BS2=BS1  ^ 

;-• . 

C$2=(rBl 

. 0S2=DBT . . . . . . . . 

ES2=EST . - . . . . - . 

c 

FS2^=r:$T 

AL2=ALI 

Bt"2=BCF 

CL2=CLI . . 

. OL2SDL1 . . . . 

c 

ErC2=Er  1 

ORSTo=s  (MF-l )  »nP 

c 

■POsT./DiTsTO  '  ^  ^ 

(50  =  0] . 

A  Sl  =  r7—QT  ■**’DRB  A 

BS1=1 .*Q0*2.»DRS6 

CSl=nRS6*P1 

DSl  sOR-? ,  <J'nP‘56'»P0  . .  "  -  - - 

ESl=nPS6 


•A,. 


F■T^*S.5 


lT/12/flO 
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C 

ALl*! .♦DP?*P1+CS1*Q1*DH2/ASI 

. . -  ■  PL  1  =l  . -0R2*P'0  +DS1  «DR2»Q r/A51 — 

CLl=0R2«O0+PSl*DR2*qi/ASl 
DL1sDR2  +  ESI«'DR'C*01/AS1 
eO'=dR2  +  FSl>'DR2*Ql  /ASl 

c 

GAmmas-A«;1<»al2/CL2/ALPHA  (MF  +  2) 

Cs- A  s  i  »  Ar2/CL2  .  ■“ 

AsAS1.»BL?/CL2+RS1*AL1/CL1-CS1 

NlsNF*l 


-AlTPHAl^n ')  VA  -G  A  M  WA 
B(Nl)=DSl-BSI»flLl/CLl 
SO=FSl-0Fl»eLl/CLI 


-T['^  S 1  *  E  L  ?  /  cL?VH:sr=9S-I^15CI7CLl  - - - - - —————————— 

S2=ASI*DL?/CL2 

P2*S2-*eP  (  WF+2tITV)  ♦$!*£«  (NF  +  l  f  lTV)  *S0*ER(NF»  ITV) 

st(i^’i)=e?-gavima*si  rNi-nr 

. OPGl  (Ml'+  j  ')  =AL2/CL2 . ; . . . .  . . . . 

ORG2(Nl  +  n=-BL's:/CL2 

ORGT=:DL?»ER  (Ml  t-l  t  ITV)  /CL^+XLZ^ER  (Nl  ♦  ITV)  /CU?  ^ 

^  .  C  .  .  .  ^  L./' 

2010  C0^4TI^^UE 

2:; 

C 

C  00  FTPS!  POIMT 

-  ,  .  c 

■  ■  .  PlsPOf:  -V 

. 

Ql=Qn 

Bs  j=B  Sl'  . . . . . . . . . . . . . . . 

j::  : 

CS2SC51  T—  ^  ^  ^  . 

ES2*ES1  , 

,  ■  FS2sFSl,^ 

Ai:2=TL  t 

. BL2*PL1 . . 

CL2=CLl  . ^ . . 

DL2=0L1 

EL2SEL1 

P0  =  0,. 

00=^1/?, 

ASl  sl  .-OTWDPFG . . .  . 

- 

eM=l  ,+Qn*2,'"UP56 

CSlsDRS6»Pl 

DSl *DP-2.»0RS6«P0 

tM=nR55 

. FSl='2*»no$(S . . . - . . 

. AL  T-r.  +nD?*oi'  *CS1»Q1»DH2/AS1 . 

4 

HLI  =  1  .-nP2«pn-+05i*DRH»^QT/A5I  ----- 

CLl =nP2»00+BSl»DR2*01/ASl 

DLlsnP2+FFl»DP£*Ql/ASl 

t:Li=np?^FSl^nR;?<^Qi/A5i  ”  ' 

C 
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C5:-A?;i»Al?/CL2 
GflMMA=-ASl«AL?/CL2/ALPHA (3) 
A=ASI»PL?/CL?*PS1«AL1/CL1-CS1 

- A-cpmrr?  )  s  A -G  A  mm  a  *  b  f  3  r - ^ - 

B(2)=y. 

S0=FSl-e5l*BLl/CLl 

s  1  =  A  s  1 «  E  i:?/  c  L  2  ♦  ES I  -8  ST^OUT  /^CUI 

S2=AS1*DL2/CL? 

F?sS?»ER(3»ITV)  +S1*ER  (ZtlTV)  ♦SO^^ER  ( 1  f  ITV) 

■  0P6  r  ( 3T = AL2/ci:2  ~ 

0RG2(3)s-BL^/CL2 

OFG3 (3) 5DL?«ER (3* ITV) /CL2*Et2*ER (2» ITV) /CL2 
OPG I  ("3 )  =Ta L  r/CL  1 "  "  “ 

0PG2 (2) s-BLA/CLl 

0PG3{?)=Dl.l»ER(2»ITV)/CLl+ELl*ER(2trTV)/CLl 
STC2T^PZ-ryKmii*SU  31  ~ 

ER(2»rTV)=sr (2)/ALPHA <2) 

c  ____  _  _  _ .  _ _ _ _ _  _ _ _ _ _ 

c 

C  _ .  . 

c 

.  DO  2020  Is3»MR 
2020  CONTTMI.IE'* . 

ERClflTVjsO.  . 


rp"mfm  TTVT^.~ 
P»^I  (MPlf  TTV)  =0  , 


00  POP?  r=i »Mo 


- - (VJF-NJPI  .j - ^ - — - - — - - 

PHI  (WF»ITV)=0PS1  eNF*l)*ER(NF-*.l,ITV)  ♦ORG2(NF*l)*EH(NFtrTV)+0RG3( 
iiNF+l  ) 

- ezcm  p'^ttv  r=E\r!n^rrMF  *  i  t  vt 

2022 CONTINUE . 

pftijrm . . 

■■  END 


Appendix  IV 

Abstract  of  paper  on  experimental  results  given 
at  the  33rd  Gaseous  Electronics  Conference. 
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Norman,  Oklahoma 

DO  NOT  WRITE  IN  THIS  SPACE 

AUTHORS  PLEASE  NOTE 

E.E.  Kunhardt 

No 

Indicate  in  what  subject  category  your 
paper  best  belongs  (see  preliminary 

Ionized  Gas  Laboratory 

Arrpptpri!  Yp5  No  _ 

Texas  Tech  University 

Session  _  _ 

announcement  for  category) 

P.O.  Box  4439 

Numher 

1  2  3  4  5  6 

©  8  9  10  11  12 

Lubbock,  Texas  79409 

Date  Conf-  _ 

13  14  15 

Nanosecond-Pulse  Breakdown  in  Gases  at  High  Over¬ 
voltages  -  E.E.  KUNHARDT,  S.  LEVINSON,  and  M.  ALLEY, 

Texas  Tech  University*  —  Physical  insight  into  the 
breakdown  initiation  processes  and  the  formation  time  of 
nanosecond-pulse  discharges  in  a  homogeneous-field  gap 
may  be  inferred  from  the  study  of  the  observational  de¬ 
lay  time.  In  this  paper,  we  report  the  results  of  ex¬ 
periments  carried  out  to  measure  the  observational  delay 
time  as  a  function  of  gas  pressure  (up  to  1600  Torr)  and 
percent  overvoltage  (up  to  1200%).  The  experiments  were 
carried  out  in  an  evacuable  spark  gap  chamber  that  forms 
part  of  a  coaxial  line  system.  Trapezbiday  voltage 
pulses  in  the  range  of  50  to  160  kV  and  50  ns  duration, 
with  4  ns  risetime  were  applied  to  the  gap.  The  elec¬ 
trode  separation  can  be  varied  up  to  3  cm.  Briefly,  the 
results  indicate  that  there  is  a  minimum  in  the  mean  de¬ 
lay  time  as  a  function  of  percent  overvoltage.  This  im¬ 
plies  that  there  are  conditions  of  overvoltage  for  which 
avalanche  development  is  optimum.  The  implication  of 
these  results  as  to  the  nature  of  the  pre-breakdown  pro¬ 
cesses  will  be  discussed. 

*  This  work  was  supported  by  the  Naval  Surface  Weapons 
Center,  Dahlgren,  Va.  under  Contract  No.  N60921-79-C- 
A187. 


61 


References 

1.  Troian,  J.  K. ,  Charbonnier,  F.  M.,  Collins,  F.  M.,  and  Guenther,  A.  H., 
Exploding  Wires  III,  pp.  361-389  0964). 

2.  Dakin,  T.  W.,  Luxa,  G.,  Opperman,  G.,  Vigreux,  J.,  Wind,  G.,  and 
Winkelnkemper,  H.,  Electra  61  (1974). 


